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ABSTRACT 

This  paper  deals  with  the  applicability  of  empirical  para¬ 
meterizations  of  stratospheric  transports  of  chemically  inert 
tracers  to  predictive  or  prognostic  two-  and  one-dimensional 
photochemical  models  of  the  stratosphere  and  troposphere  for 
the  forecasting  of  anthropogenic  effects  on  atmospheric  ozone. 

The  scope  of  this  paper  therefore  includes  (1)  a  critical  re¬ 
view  and  assessment  of  the  prognostic  utility  of  the  parent 
Reed  and  German  (1965)  2-D  parameterization  of  stratospheric 
transports,  (2)  the  implied  assumption  in  representative  sub¬ 
sidiary  2-D  parameterizations  used  or  for  use  in  2-D  photochemical 
models,  (3)  use  of  GCM/tracer  model  data  for  a  chemically  inert 
tracer  for  the  assessment  of  the  prognostic  utility  of  both  2-D 
and  1-D  parameterizations  of  stratospheric  transports,  and  (4) 
the  outlook  for  the  development  of  prognostic  parameterizations 
of  stratospheric  transports. 
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SUMMARY 


This  paper  deals  with  the  applicability  of  empirical  para- 
meterizations  of  stratospheric  transports  of  chemically  inert 
tracers  to  predictive  or  prognostic  two-  and  one-dimensional 
photochemical  models  of  the  stratosphere  and  troposphere  for 
the  forecasting  of  anthropogenic  effects  on  atmospheric  ozone. 
Because  the  prognostic  value  of  available  two-dimensional  (2-D) 
parameterizations  of  stratospheric  tracer  transports  depend  on 
the  viability  of  the  parent  Reed  and  German  (1965)  formulation, 
emphasis  is  placed  here  on  basic  limitations  of  such  formulation. 
The  prognostic  viability  of  1-D  parameterizations  of  strato¬ 
spheric  tracer  transports  is  evaluated  from  comprehensive  3-D 
GCM#/tracer  model  data  for  a  chemically  inert  tracer  as  a  func¬ 
tion  of  tracer  configuration  in  both  the  space  and  time  domains 
as  well  as  of  tracer-source  location. 

General  and  specific  results  for  both  the  parent  and  sub¬ 
sidiary  2-D  parameterizations  of  stratospheric  tracer  transports 
are  as  follows: 

•  The  2-D  Reed  and  German  parameterization  of  tracer  trans¬ 
ports  in  the  critical  region  of  the  lower  stratosphere 
is  based  on  the  crucial  assumption  that  the  macroscale 
vertical  and  meridional  eddy  transports,  driven  by  well- 
structured  stratospheric  wave  motions,  may  be  described 
by  the  same  turbulent  mixing  mechanisms  as  in  the  tropo¬ 
sphere.  As  a  consequence  of  this  assumption,  there  is 
a  need  for  drastic  physical  reinterpretation  of  the 
parent  Reed  and  German  formulation  of  stratospheric  eddy 
transports,  a  formulation  that  is  based  on  an  adaptation 

1 - 

General  circulation  model. 
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of  the  Prandtl  mixing-length  hypothesis  for  microscale 
turbulence.  Hence,  the  applicability  of  the  subsidiary 
2-D  parameterizations  in  available  photochemical  models 
for  the  forecasting  of  anthropogenic  ozone  effects  is 
open  to  question. 

•  Viable  prognostic  parameterizations  of  stratospheric 
eddy  transports  for  2-D  photochemical  models  require, 
in  addition  to  the  foregoing  reinterpretation  of  the 
physical  nature  of  the  eddy  fluxes,  drastic  modifica¬ 
tions  of  the  parent  Reed  and  German  formulation  in 
regard  to  (1)  the  calculation  of  the  statistics  of  the 
zonal  average  of  both  the  slopes  of  the  mixing-length 
paths  or  surfaces  (a)  and  their  deviation  from  such 
average  (a*2),  and  (2)  the  basic  characteristics  of  the 
eddy  transport  coefficients  (K's)  which  (a)  may  not 
have  the  Kyz  =  symmetry  as  derived  from  implicit 
assumptions  in  the  parent  formulation  (Eqs.  7  to  10  vs 
Eqs.  15  to  19)  and  (b)  are  not  independent  of  the  zonal 
(u)  and  meridional  (w,  v)  circulations  as  assumed  in 
such  parent  formulation. 

•  A  strict  interpretation  of  the  Reed  and  German  2-D 
formulation  of  eddy  transports  leads  to:  (1)  apparent 
undeterminability  from  vertical  and  meridional  GCM  wind 
data  (Tables  5  through  7)  of  the  basic  a*2/ci2  statistics 
of  the  slope  (a)  of  mixing-length  paths  and  (2)  over¬ 
specification  of  two  eddy  transport  equations  (Eqs.  4,  5) 
for  one  unknown  eddy  transport  coefficient  (K  )  as  a 
consequence  of  the  required  equality  of  the  Kyz  and  K„y 
eddy  transport  coefficients  (Eq.  9).  This  latter  result 
is  opposite  to  current  practice  of  using  an  undetermined 
system  of  two  eddy  flux  equations  (4,  5)  for  determining 

three  unknown  eddy  coefficients  (K  ,  K  ,  and  K  :  with 

yy  yz  zz’ 
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•  Since  the  2-D  parameterizations  of  tracer  transports 

used  or  for  use  in  current  representative  photochemical 

models  (Tables  1  to  3)  assume  equality  of  the  K  „  and 

yz 

K  coefficients  in  the  stratosphere  and  ignore  the  de- 
zy 

pendence  of  such  coefficients  on  the  zonal  and  meridional 
circulations,  they  cannot  therefore  rely  on  the  viability 
of  the  Reed  and  German  formulation  to  justify  their 
prognostic  utilization.  Moreover,  it  is  found  that 
they  are  based  on  implied  assumptions  of  the  basic 
a*2/a2  statistics  that  are  (1)  not  always  positive,  even 
though  they  involve  a  ratio  of  quadratic  terms  (Table  4); 
(2)  drastically  different  among  themselves  (Table  4) 
instead  of  remaining  fixed  by  the  (same)  wind  variability 
of  the  stratosphere  (Eq.  6);  and  (3)  inconsistent  with 
the  parent  Reed  and  German  formulation  due  to  the  apparent 
undeterminability  (Eq.  14)  of  such  statistics. 

•  Use  of  3-D  GCM/tracer  model  data  for  the  evaluation  of 
2-D  parameterizations  of  eddy  transports,  when  they  are 
based  on  the  assumption  that  a*2/ci2  =  0  everywhere  in 
the  lower  stratosphere  at  middle  and  high  latitudes 
(Table  1),  indicate  the  following:  (a)  the  eddy  co¬ 
efficients  would  have  to  become  negative  in  some  regions 
of  the  stratosphere  (Tables  9  and  10),  and  (b)  when  air¬ 
craft  emissions  are  introduced  continuously  in  the  lower 

stratosphere,  the  a  (or  K  /K  )  parameter  would  also 

v  **  y  % y 

have  to  change  instead  of  remaining  fixed  by  the  wind 
variability  (Table  8). 

•  The  significant  differences  in  the  magnitude  of  the 

eddy  transport  coefficients,  as  given  by  available  repre¬ 
sentative  2-D  parameterizations  of  eddy  transports  in 
the  lower  stratosphere  (e.g.  Tables  2  vs  3),  appear  to 
stem  from  two  factors:  (1)  use  of  different  types  of 
observed  tracer  data  for  the  determination  of  such  co¬ 
efficients  and  (2)  the  fact  that,  as  indicated  by 
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GCM/tracer  model  data,  the  eddy  transports  of  anthropogenic 
tracers  depend  on  the  time  (age  in  years)  and/or  space 
(source  location)  configurations  of  the  tracers  in  tne 
stratosphere . 

•  Available  preliminary  results  from  a  Lagrangian  instead 

of  Eulerian  description  of  air  parcel  motions  caused 

by  a  planetary  wave  in  the  stratosphere  indicate  the 

possibility  for  a  prognostic  2-D  formulation  of  eddy 

transports  which  would  (1)  have  no  need  for  the  Reed 

and  German  constraint  on  an  equality  of  the  Kyz  and 

K  eddy  transport  coefficients,  and  (2)  relate  the 
zy 

eddy  transport  coefficients  to  the  period  of  oscillation 
of  air  parcels  due  to  planetary  waves,  which,  in  turn, 
relates  to  the  zonal  circulation. 


In  regard  to  the  prognostic  applicability  of  empirical  1-D 
parameterizations  of  stratospheric  transports  in  1-D  photochemical 
models,  comprehensive  3-D  GCM/tracer  model  data  (Table  11)  indi¬ 
cate  the  following: 


•  The  magnitude  of  the  1-D  transport  coefficients  in  the 
critical  region  of  the  lower  stratosphere  can  (1)  be¬ 
come  negative  for  SST  emissions,  a  result  that  implies 
distortions  of  the  K(z)  profile  and  (2)  depend  on  the 
tracer  configuration  in  the  stratosphere  as  given  by 
the  source  location,  source  characteristics,  and  tracer 
age  or  years.  Therefore,  the  empirical  parameteriza¬ 
tions  of  stratospheric  transports  used  in  1-D  photochemical 
models  appear  to  lack  physical  basis  for  their  prognostic 
utilization,  since  they  are  assumed  to  be  always  positive 
and  independent  of  tracer  configuration  in  the  strato¬ 
sphere. 


Although  the  magnitudes  of  1-D  coefficients  derived  from 
3-D  GCM/tracer  model  data  may  be  dependent  on  the  vertical  re¬ 
solution  of  the  model  in  the  critical  tropopause  and  lower 


stratosphere  regions,  the  trends  established  by  extensive  cal¬ 
culations  of  such  a  model  during  the  last  several  years  may  not 
be  ignored  in  the  absence  of  worldwide  observations  of  suitable 
tracers  for  an  ensemble  of  many  years.  In  fact,  a  promising 
avenue  to  derive  more  exact  1-D  transport  coefficients  is  to 
improve  the  vertical  resolution  of  the  GCM/tracer  model  instead 
of  ever  attempting  to  collect  empirical  data  for  a  suitable 
tracer  in  the  lower  stratosphere  over  a  worldwide  space  domain 
and  for  an  ensemble  of  many  years. 
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I.  INTRODUCTION 


A  need  to  forecast  long-term  effects  of  high-altitude 
aircraft  engine  emissions  (N0x  and  HgO)  on  atmospheric  ozone 
as  a  function  of  altitude,  latitude,  and  season  has  led  to  an 
interest  in  the  development  of  predictive  or  prognostic  two- 
dimensional  (2-D)  photochemical  models  of  the  stratosphere  and 
troposphere  (e.g.  Hidalgo  and  Crutzen,  1977;  and  Widhopf,  Glatt, 
and  Kramer,  1977).  A  basic  component  of  such  models  is  the 
description  of  large-scale  vertical  and  latitudinal  or  meridional 
transports  of  tracer  mixing  ratio  for  the  determination  of  the  2-D 
distributions  of  relevant  natural  and  anthropogenic  chemical 
species.  The  development  of  prognostic  2-D  photochemical  models 
of  the  stratosphere  and  troposphere  has  then  consisted  in  the 
generation  of  numerical  solutions  of  time-dependent  partial 
differential  equations  for  the  conservation  of  each  of  the 
involved  chemical  species  with  specified  (Hidalgo,  1978):  (1) 

vertical  and  meridional  mean  winds  as  well  as  eddy  transport 
coefficients  as  a  function  of  altitude,  latitude,  and  month 
(season);  (2)  photolysis  rates  and  short-wave  solar  radiation 
for  the  relevant  photochemistry;  (3)  reaction  rates  for  the  in¬ 
volved  chemical  kinetics;  (4)  rate  constants  for  the  heterogeneous 
chemistry  involving  sinks  of  water-soluble  chemical  species  in 
the  lower  troposphere;  and  (5)  anthropogenic  sources  of  chemical 
tracers  as  a  function  of  atmospheric  location  and  time  (years). 

The  computing  scheme  is  completed  by  using  for  each  of  the 
chemical  species  initial  conditions  throughout  the  model  atmos- 
sphere  as  well  as  appropriate  boundary  conditions  at  the  earth's 
surface  and  stratopause. 
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Two  main  applications  of  2-D  photochemical  models  during 
the  last  several  years  have  been  (1)  diagnostic  calculations 
of  the  prevailing  mixing  ratio  of  relevant  chemical  species  as 
a  function  of  altitude,  latitude,  and  season  in  the  stratosphere 
and  troposphere;  and  (2)  prognostic  forecasts  of  long-term, 
global-scale  ozone  effects  from  specified  future  aircraft  emis¬ 
sions  as  a  function  of  altitude  and  latitude.  Since  effective 
regulations  designed  to  prevent  predicted  long-term  adverse 
ozone  effects  on  a  global  scale  would  require  international 
participation,  the  scrutiny  and  documentation  of  each  of  the 
above  foundations  of  2-D  photochemical  models  must  be  such  as 
to  promote  international  credibility  of  their  predictions  con¬ 
cerning  any  anthropogenic  threat  to  atmospheric  ozone. 

The  objective  of  this  paper  is  to  examine  the  description 
of  the  large-scale  vertical  and  meridional  transports  of  chem¬ 
ical  species  as  used  or  intended  for  use  in  current  representative 
2-D  photochemical  models  of  the  stratosphere  and  troposphere. 
Furthermore,  since  such  descriptions  of  tracer  transports  are 
based  on  extensions  or  interpretations  of  the  parent  Reed  and 
German  (1965)  formulation  for  chemically  inert  tracers  in  the 
stratosphere,  emphasis  must  be  placed  here  on  this  parent  formu¬ 
lation.  Considerations  of  interest  are  then  on  the  following: 

•  Need  to  distinguish  between  the  diagnostic  use  of 
empirical  parameterizations  of  stratospheric  trans¬ 
ports  and  the  prognostic  use  of  such  empiricism  for 
the  forecasting  of  future  anthropogenic  effects  on 
ozone . 

•  Need  to  compare  on  a  common  basis  different,  repre¬ 
sentative  2-D  parameterizations  of  tracer  transports 
used  or  for  use  in  2-D  photochemical  models. 

•  Need  to  bring  out  any  basic  limitation  in  the  pioneer¬ 
ing  Reed  and  German  formulation,  which  was  devised 
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over  a  decade  ago  for  diagnostic  calculations  of 
chemically  inert  tracers  in  the  stratosphere  and 
in  the  absence  of  current  3-D  general  circulation 
model  (GCM)  data  for  unmeasurably  small  large-scale 
vertical  transports  in  the  stratosphere.  An  indirect 
tribute  to  the  difficulty  of  such  a  pioneering 
effort  is  the  fact  that  there  is  not  yet  an  adequate 
solution  to  the  2-D  formulation  of  tracer  transports 
in  the  stratosphere. 

•  Need  to  clarify  past  interpretations  concerning 

basic  concepts  of  the  parent  Reed  and  German  formu¬ 
lation  of  tracer  transports  in  the  stratosphere. 

Hence,  the  main  scope  of  this  paper  includes  the  following: 

(1)  review  of  the  Reed  and  German  parameterization  of  the  transports 
of  chemically  inert  tracers  in  the  stratosphere;  (2)  implied 
assumptions  in  the  parameterizations  of  tracer  transports  in 
current  representative  2-D  photochemical  models;  (3)  apparent 
undeterminability  of  fundamental  statistics  in  the  Reed  and  German 
parameterization;  (4)  constraints  in  basic  Reed  and  German  results 
that  are  carried  over  to  the  2-D  parameterizations  in  photo¬ 
chemical  models;  (5)  past  interpretations  of  basic  Reed  and 
German  results;  and  (6)  outlook  for  future  progress  in  the  de¬ 
velopment  of  prognostic  2-D  formulations  of  tracer  transports. 
Finally,  because  of  the  widespread  use  of  1-D  parameterizations 
of  atmospheric  transports  in  the  more  versatile  1-D  photochemical 
models,  this  paper  includes  considerations  of  the  limitations 
of  1-D  empiricism  for  the  prognostic  representation  of  atmospheric 
transports. 
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II.  REVIEW  OF  REED  AND  GERMAN  FORMULATION 

The  basic  aim  of  2-D  parameterizations  of  stratospheric 
transports  in  photochemical  models  is  to  describe  the  zonal- 
monthly  (seasonal)  statistical  averages  of  vertical  and  meridional 
transports  of  chemical  species  in  stratospheres  that  are  character¬ 
ized  by  either  the  absence  or  presence  of  aircraft  engine 
emissions  of  N0x  (NO,  NO^)  and  H^O.  Because  of  the  interest  in 
monthly  (seasonal)  anthropogenic  ozone  effects,  the  primary 
concern  is  with  the  statistics  of  tracer  transports  by  large- 
scale  wave  motions.  Two  basic  assumptions  used  ;.n  representative 
2-D  photochemical  models  of  the  stratosphere  and  troposphere  are 
(1)  decoupling  of  the  chemistry  from  the  temperature  and  wind 
fields,  and  (2)  use  of  the  Heed  and  German  adaptation  of  the 
Prandtl  (1925)  mixing-length  hypothesis  to  the  macroscale  eddy 
transports  in  the  stratosphere.  Each  of  these  assumptions  is 
described  below. 

A.  DECOUPLING  OF  CHEMISTRY  FROM  THE  TEMPERATURE  AND  WIND  FIELDS 

The  first  basic  assumption  in  representative  2-D  photochem¬ 
ical  models  of  the  stratosphere  and  troposphere  is  that  the 
conservation  equations  for  chemical  species  are  decoupled  from 
the  primitive  equations  of  the  general  circulation  (e.g.  Lorenz, 
1967).  Since  the  latter  equations  can  be  used  to  obtain  esti¬ 
mators  of  the  monthly  (seasonal)  statistics  of  the  temperature 
and  wind  fields  in  the  stratosphere,  an  implicit  assumption  in 
the  2-D  photochemical  models  is  that  the  anthropogenic  redis¬ 
tributions  of  radiatively  and  chemically  active  tracers  (e.g. 
ozone)  will  not  produce  first-order  effects  on  the  monthly 
statistics  of  the  temperature  and  winds  in  the  stratosphere.  A 
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consequence  of  such  decoupling  is  that  such  2-D  photochemical 
models  cannot  then  take  into  account  feedback  effects  of  the 
calculated  ozone  redistributions  on  the  original  statistics  of 
the  background  temperature  and  wind  fields.  Obviously,  the 
validity  of  the  decoupling  of  the  chemistry  from  the  temperature 
and  wind  fields  would  become  open  to  question  as  the  magnitude 
of  the  calculated  changes  in  ozone  by  all  anthropogenic  effects 
become  increasingly  large. 

The  Interactions  of  the  photochemistry  and  chemical  kinetics 
with  the  temperature  and  wind  fields  are  given  by  a  set  of  con¬ 
tinuity  equations  for  the  chemical  species  in  an  oxygen-hydrogen- 
nitrogen-chlorine  atmosphere,  i.e.,  by 


(1) 


where  R^  is  the  mixing  ratio  for  any  ith  of  N  given  chemical 
species,  t  the  time  and  P.^  as  well  as  Li  the  corresponding  pro¬ 
duction  (including  anthropogenic  sources)  and  loss  (or  sink) 
chemical  mechanisms.*  The  left  side  of  Eq.  (1)  contains  the 
influence  of  motions  as  identifiable  by  expanding  the  total 
derivative  dR^/dt. 

The  2-D  continuity  equations  for  the  conservation  of  chemical 
species  are  obtained  by  (1)  taking  the  zonal-monthly  statistical 
average  of  Eq.  (1)  for  an  ensemble  of  many  years  and  (2)  splitting 
the  transports  of  mixing  ratio  into  those  by  the  mean  winds  and 
those  by  wave  or  eddy  motions.  For  the  case  of  chemically  inert 
tracers,  as  in  the  Reed  and  German  formulation,  =  0, 

=  R  and  Eq .  (1)  simplifies  to  dR/dt  =0.  In  spherical 

« - 

The  terms  P^  and  Lj_  in  Eq .  (1)  have  absorbed  the  air  number 
density  n  as  obtained  from  use  of  the  continuity  equation  for 
both  the  ith  chemical  species  (ni)  and  air  (n),  together  with 
Ri  =  nj/n. 


6 


coordinates,  the  zonal-monthly  statistical  average  of  dR/dt 
becomes  as  in  the  Reed  and  German  Eq.  (2),  i.e., 

||  =  -  V-VR  -  i  |j  o  OTT  -  p  OPS9  ,  (2) 

where  V  is  the  2-D  wind  vector,  V  the  corresponding  gradient 
vector  for  the  mixing  ratio,  p  the  air  density,  w  and  v  the 
upward  and  poleward  winds,  a  is  the  mean  radius  of  the  earth, 
z  the  altitude,  and  0  the  latitude.  The  Reed  and  German  notation 
(")  denotes  the  estimators  for  the  zonal-monthly  ensemble  aver¬ 
age  of  (  ) ,  whereas  (  )*  is  the  deviation  of  (  )  from  such 
average;  i.e.  (  )#  =  (  )  -  (“).  The  right-hand  side  of  Eq.  (2) 
exhibits  the  roles  of  the  transports  of  mixing  ratio  by  (1) 
the  meridional  circulation  (V*7R)  or  average  vertical  and 
meridional  winds,  (2)  the  vertical  eddy  motions  (2nd  term)  and 
(3)  the  meridional  eddy  motions  (remaining  term). 

B.  ADAPTATION  OF  MIXING  LENGTH  HYPOTHESIS  TO  THE  STRATOSPHERE 

The  numerical  solution  of  the  continuity  equation  (2)  for 
a  chemically  inert  tracer  requires  elimination  of  the  eddy 
transport  terms  through  their  correlations  with  the  vertical 
and  meridional  components  of  the  gradient  vector  VR.  For  this 
purpose,  the  2-D  Reed  and  German  representation  of  eddy  trans¬ 
ports  adapts  the  Prandtl  mixing-length  hypothesis  for  turbulent 
transports  of  momentum  in  microscale  flow  to  the  macroscale 
transport  by  wave  motions  of  chemically  inert  tracers  in  the 
stratosphere.  As  a  consequence  of  this  assumption,  the  Reed 
and  German  formulation  carries  over  to  the  stratosphere  the  im¬ 
plicit  assumptions  that  (a)  the  eddy  tracer  transports  in  Eq.  (2) 
are  independent  of  both  the  zonal  and  meridional  circulations, 
when  in  fact  they  are  intimately  related  (e.g.  Hunt  and  Manabe, 
1968;  Holton,  1975);  and  (b)  the  eddy  transports  can  then  be 
correlated  with  gradients  of  mixing  ratio,  as  in  the  Prandtl 
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mixing  length  hypothesis.  Thus,  if  x  denotes  a  dependent  var¬ 
iable  such  as  potential  temperature  or  ozone,  a  direct  application 
of  the  Prandtl  mixing  length  hypothesis  would  give  v#x*  * 

-  Ky  3x/3y;  where  y  =  a0  is  taken  as  positive  (Ay  >  0)  in  the 
poleward  direction.  However,  a  basic  difficulty  then  was  that 
for  a  poleward  eddy  flux  (i.e.  v*x*  >  0),  with  K  >  0,  then 
3x/9y  <  0J  i.e.,  the  poleward  eddy  flux  had  to  take  place  always 
in  a  downgradient  direction  from  high  to  low  values  of  x  in 
opposition  to  observed  countergradient  eddy  fluxes  of,  for  ex¬ 
ample,  ozone  ( v*03# )  or  heat  . ( v#T# ) .  Hence,  in  order  to  adapt 
the  Prandtl  mixing  length  hypothesis  to  countergradient  eddy 
fluxes  from  low  to  high  values  of  x,  it  was  postulated  that  the 
slope  (a)  of  mixing  length  paths  or  surfaces  would  have  to  be 
larger  than  the  slope  (6)  of  the  tracer  isopleths  x  or  isentropes 
when  x  represents  the  potential  temperature.  This  concept  is 
illustrated  in  Fig.  1,  as  given  by  Reed  and  German.  The  basic 
results  of  this  2-D  formulation  of  the  eddy  transports  are 
then  as  follows: 


X*  -  -  ■  -  «y  &  +  H>  t3> 

where  t  is  the  mixing  length  vector  (l  -  100  km)  with  horizon¬ 
tal  and  vertical  components  l „  and  £_.  When  x  represents  mixing 

y  z 

ratio  in  Eq.  (3),  the  desired  correlations  of  the  horizontal 
and  vertical  eddy  fluxes  in  Eq.  (2)  with  the  components  of  VR 
are  given  by 


v*R* 

w*R* 

where  Kyy  I  7*^,  Kyz 
the  fact  that  a  «  10~H 
length  direction: 


3R  ,  j.  3Rv 

■yy  3y  yz  * 

(4) 

3R  i  v  3S \ 

zy  "57  zz  "Tz  * 

(5) 

Kzy  =  *nd  Kzz  - 

Using 

*  -  (K 

-= 

and  that  the  wind  V  is  along  the  mixing 
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FIGURE  1.  Schematic  representation  of  the  slope  (a)  of  the  mixing  length 
path  necessary  for  countergradient  eddy  fluxes  in  the  Prandtl 
mixing  length  hypothesis,  as  given  by  Reed  and  German  (1965). 

The  eddy  flux  is  countergradient  here  because  both  3x/3y  and 
v*x*  are  positive.  The  former  condition  takes  place  because 
Ax  >  0  when  Ay  >>  0,  and  the  latter  because  the  eddy  flux  is 
poleward;  a  condition  that  can  be  visualized  by  letting  x  repre¬ 
sent  the  potential  temperature.  The  dashed  arrows  then  indicate 
the  mixing  at  the  origin  of  northward  moving  air  parcels  that 
are  warmer  than  ambient  with  southward  moving  air  parcels  that 
are  colder  than  ambient.  Note  that  the  countergradient  eddy 
flux  of  heat  v*T*  is  poleward  whenever  a  >  6.  Similar  argu¬ 
ments  apply  when  x  represents  ozone  to  yield  a  countergradient 
flux  characterized  by  3$3/3y  >  0  and  v*0j*  >  0. 
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v*  =  V  cos  a  -  V 
w#  =  V  sin  a  -  Va 


(6) 


=  £  cos  a  -  l 
l  -l  sin  a  =  la 

With  Eqs.  (6),  the  eddy  coefficients  K  .  K  .  K  .  and  K  can 

yy  yz*  zy  zz 

be  expressed  as  follows: 


Kyy 

=  vF 

eosza  ®  VI  »  K 

(7) 

Kyz 

=  vi 

sinacosa  -  Via 

«  a  K 

(8) 

Kzy 

■  ^ 

=  vF 

sinacosa  «  Via 

Kyz 

(9) 

Kzz 

=  w  *1 
z 

=  vF 

sin2a  -  Via 2  ~ 

(a2  +  a*2)  K 

(10) 

where  it  is  aesumed  that  a  and  a*  are  independent  of  V  and  l, 
i.e.,  that  they  do  not  depend  on  the  wind  speed  or  the  length 
of  the  mixing  path.  Substitution  of  Eqs.  (7)  through  (10)  into 
(M  and  (5),  using  6  =  tan  B  =  -  Ol/9y )/OE/3z ) ,  alternative 
expressions  for  the  eddy  fluxes  are  given  by 


=  -  Kyy  (1  -  !}  § 


(id 


w*R* 


Kzz  (1 


SB  v  3R 
a2"+"a«4; 


(12) 


Equation  (11)  indicates  that  (a)  when  a/B  =  0,  the  horizontal 
eddy  flux  becomes  downgradient  as  in  the  Prandtl  formulation  for 
the  microscale  eddy  transport  of  momentum;  (b)  when  a/B  =  1,  the 
horizontal  eddy  flux  vanishes;  and  (c)  when  S/B  >  1,  the  hori¬ 
zontal  flux  becomes  countergradient.  Similar  consideration  can 
be  made  to  the  vertical  eddy  flux  as  given  by  Eq.  (12). 
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III.  PARAMETERIZATION  OF  EDDY  TRANSPORTS  FOR 
2-D  PHOTOCHEMICAL  MODELS 


The  continuity  equations  in  representative  2-D  photochemical 
models  are  as  Eq.  (2),  except  that  R  must  be  replaced  by  R^  and 
the  right-hand  side  must  incorporate  the  and  terms.  The 
use  of  Eqs.  (4)  and  (5)  for  each  ith  chemical  tracer  allows 
elimination  of  the  eddy  terms  in  these  equations,  which  become 
differential  equations  for  the  mixing  ratio  of  each  tracer  with 
specified  average  meridional  winds  and  eddy  transport  coeffi¬ 
cients.  A  fundament&l  approach  in  2-D  photochemical  models 
then  consists  in  (1)  assuming  that  the  eddy  transport  coeffi¬ 
cients  are  independent  of  the  chemical  tracer,  i.e.,  that  they 
depend  only  on  altitude,  latitude,  and  month  for  a  specified 
meridional  circulation  and  (2)  establishing  the  magnitude  of 
the  eddy  transport  coefficients  from  observations  of  distributions 
of  natural  (i.e.  ozone,  water  vapor)  and/or  anthropogenic  (i.e. 
radioactive  debris  from  past  nuclear  explosions)  tracers  in 
the  lower  stratosphere.  The  parameterizations  of  eddy  transports 
have  then  been  reduced  to  the  use  of  two  eddy  tracer  flux 
equations  (4,  5)  for  three  unknown  eddy  transport  coefficients 
(Kyy,  Kyz,  Kzz,  i.e.  with  Kzy  =  Kyz)  at  given  stratospheric 
meridional  locations  and  month  (season).  This  procedure  has 
therefore  required  the  use  of  additional  constraints  so  as  to 
resolve  the  undetermined  foregoing  system. 

Three  representative  2-D  parameterizations  of  eddy  trans¬ 
ports  utilizing  nearly  the  same  meridional  circulation  (i.e. 

Louis,  1974)  are:  (1)  that  in  the  Crutzen  photochemical  model 
(Hidalgo  and  Crutzen,  1977),  (2)  that  in  the  Widhopf,  Glatt, 
and  Kramer  (1977)  model  and  (3)  that  given  by  Danielsen  and 
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Louis  (1975).  The  Crutzen  parameterization  of  eddy  transports 
utilized  Eqs.  (7)  through  (10),  together  with  continuity  equations 
for  ozone  and  water  in  the  presssure  instead  of  the  altitude  co¬ 
ordinate  as  it  is  done  in  the  other  two  parameterizations .  The 
Wldhopf  et  al.  parameterization  utilized  Eqs.  ( 4 ) ,  (5),  and  (9), 
together  with  the  continuity  Eq.  (2)  for  chemically  inert  tracers 
and  (a)  extrapolations  to  the  stratosphere  (by  Luther,  1973)  of 
the  Oort  and  Rasmusson  (1971)  data  for  eddy  transports  of  heat 
in  the  upper  troposphere  at  latitudes  9  _<  50°  N  and  (b)  use  of 
rather  crude  observations  of  the  meridional  distributions  of  mix¬ 
ing  ratio  for  anthropogenic  chemically  inert  radioactive  tracers. 
The  Danielsen  and  Louis  parameterizations  utilized  Eqs.  (4),  (5), 
and  (9),  together  with  Eq.  (2)  for  ozone  in  the  lower  strato¬ 
sphere  on  the  assumption  that  ozone  could  be  considered  there  to 
be  a  chemically  inert  tracer;  however,  instead  of  using  the 
continuity  equation  for  water  vapor  in  the  stratosphere,  as  in 
the  Crutzen  parameterization,  the  Danielsen  and  Louis  approach 
had  to  use  constraints  that  would  insure  tropospheric  ozone 
fluxes  that  were  compatible  with  the  ozone  destruction  at  the 
ground . 

Tables  1  through  3  (Hidalgo,  1978)  illustrate  the  magnitude 
of  the  eddy  transport  coefficients  Kyy,  “  (or  Kyz/Kyy) ,  and  Kzz 
from  the  foregoing  parameterizations  at  critical  altitudes  for 
the  exchange  of  anthropogenic  tracers  between  the  lower  strato¬ 
sphere  and  upper  troposphere.  Because  of  the  significant  differ¬ 
ences  in  the  magnitudes  of  these  eddy  coefficients,  it  is  of 
interest  to  (1)  contrast  the  direct  use  of  Eq .  (10)  with  a*2  =  0 
in  the  Crutzen  parameterization  (Crutzen,  1980)  as  it  is  not  done 
in  the  other  two,  and  (2)  emphasize  the  general  dependence  of 
eddy  tracer  transports  on  both  the  space  and  time  (years)  distri¬ 
butions  of  anthropogenic  (i.e.  radioactive)  tracers  in  the 
stratosphere  (Mahlman  and  Moxim,  1978;  Hidalgo,  1978),  i.e.,  the 
use  of  the  eddy  flux  equations  (4)  and  (5)  with  Eq .  (9)  tends 
to  make  the  magnitude  of  the  eddy  coefficients  dependent  on  the 
space  and  time  (age)  distributions  of  particular  tracers  in  the 
lower  stratosphere  (Danielsen  and  Louis,  1975). 
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TABLE  1.  SAMPLE  EDDY  PARAMETERS  FOR  CRUTZEN  MODEL  DURING  SUMMER  AND  WINTER  IN  NORTHERN 

HEMISPHERE*  (Source:  Hidalgo  and  Crutzen,  T977) 
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TABLE  2.  SAMPLE  EDDY  PARAMETERS  FOR  WIDHOPF  MODELS  DURING  SUMMER  (AUGUST)  IN  NORTHERN 
HEMISPHERE  AND  WINTER  IN  SOUTHERN  HEMISPHERE*  (Source:  Widhopf,  1978) 


•The  actual  notation  for  the  eddy  coefficients  in  the  Widhopf  model  is  K^,  K  and  Kzz  (Widhopf,  Glatt  and 
Kramer,  1977).  The  angle  n.  is  given  here  by  the  ratio  K^/K^.  The  eddy  coefficients  in  this  model  are 

given  as  a  function  of  month  for  both  northern  and  southern  latitudes.  Above  data  is  applicable  for  the 
month  of  August;  i.e.,  summer  in  the  Northern  hemisphere  and  winter  in  the  Southern  Hemisphere.  See  text 
for  thje  shift  from  southern  (August)  to  northern  (February)  latitudes  for  winter. 


TABLE  3.  SAMPLE  EDDY  PARAMETERS  DERIVED  FROM  POLEWARD  TRANSPORT  OF  OZONE* 
(Source:  Danielsen  and  Louis,  1977) 
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It  should  be  noted  that  the  foregoing  use  of  the  time  depend¬ 
ent  continuity  equation(s)  for  either  chemical  or  inert  tracers 
does  not  provide  an  independent  additional  constraint  so  as  to 
fix  a  solution  for  the  three  unknown  eddy  coefficients  from  the 
two  eddy  tracer  flux  equations  (4  and  5).  As  indicated  by 
Eqs .  (11)  and  (12),  fundamental  parameters  in  the  Reed  and  German 
formulation  are  the  statistics  for  the  slope  of  mixing  length 
paths  or  surfaces  (a)  and  its  eddy  component  or  deviation  from 
this  average  (a*2).  Therefore,  the  determination  of  the  three 
eddy  coefficients  at  given  stratospheric  meridional  locations 
and  month  requires  an  additional  constraint  on  these  basic  sta¬ 
tistics  by  each  parameterization.  These  constraints  can  be 
brought  out  by  dividing  Eq .  (10)  by  the  square  of  Eq .  (8); 
after  use  of  Eq.  (7),  the  following  result  is  obtained: 


Ky.y  KZZ 
Kyz 


(13) 


Equation  (13)  provides  a  common  basis  to  compare  the  implied 
assumptions  for  the  fundamental  a* 2  statistics  in  the  foregoing 
2-D  parameterizations  of  eddy  transports.  It  is  important  to 
note  that  a*2/a2  is  controlled  only  by  the  variability  of  the 
stratospheric  winds  (Eqs.  6),  and  that  its  magnitude  must  be 
always  positive  because  it  Involves  the  ratio  of  two  quadratic 
terms.  Table  4  shows  a*2/a2  values  corresponding  to  the  three 
parameterizations  in  Tables  1  to  3.  The  results  in  Table  4 
indicate  that  these  representative  2-D  parameterizations  fail 
to'  provide  (a)  positive  values  of  a*2/a2  at  every  latitude  and 
altitude  and  (b)  nearly  identical  values  for  this  ratio  at  a 
given  location  and  season,  as  it  should  due  to  the  dependence 
of  oi*2/a2  on  the  wind  variability  of  the  same  stratosphere. 

Hence,  these  results  Indicate  that  the  foregoing  representative 
2-D  parameterizations  of  eddy  transports  are  inconsistent  both 
with  the  parent  Reed  and  German  formulation  and  among  themselves. 
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TABLE  4.  VALUES  OF  ci*2/a2  FOR  THREE  2-D  PARAMETER! ZAT I ONS  OF  EDDY  TRANSPORTS. 
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Since  the  Crutzen  model  used  a*2  *  0  (Crutzen,  1980),  then  Eq.  (13)  yields  KyyKzz/K 3Z  =  1.  The  data  in  Table  1  shows 
that  neither  of  these  two  conditions  holds  at  low  latitudes  (Hidalgo,  1978).  The  above  differences  between  the  Crut¬ 
zen  and  Widhopf  values  at  low  latitudes  may  help  explain  the  difference  in  interhemispheric  aircraft  effects  on  ozone 
from  these  models. 


A  consequence  of  the  former  result  is  that  such  parameteriza- 
tions  may  not  then  use  the  Reed  and  German  (1965)  formulation 
to  imply  physical  justification  of  their  prognostic  utility 
for  viable  forecasting  of  anthropogenic  ozone  effects. 
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IV.  APPARENT  UNDETERMINAB ILITY  OF  BASIC  a*2  STATISTICS 


The  foregoing  2-D  parameterizations  of  eddy  transports  had 

to  contend  with  the  solution  of  an  undetermined  system  of  two 

equations  (4  and  5)  and  three  unknowns  (K  ,  K  ,  K  ).  The 

yy  yz5  zz 

solution  of  such  a  system  required  then  additional  constraints 
involving  implied  assumptions  for  a*2/a2  as  given  by  Eq.  (13). 
However,  a  straightforward  alternative  procedure  would  consist 
in  a  direct  determination  of  the  a  and  a*2  statistics,  which 
would  be  derivable  from  the  Reed  and  German  Eqs.  (6),  i.e. 


tan  a 


(14) 


where  a  (0,  A,  z,  t)  is  the  slope  of  mixing  length  paths  or 
surfaces  as  a  function  of  latitude  (6),  longitude  (A),  altitude 
(z),  and  time  (t,  in  hours).  Equation  (14)  indicates  the 
following:  (1)  a  is  a  function  of  the  variability  of  the  verti¬ 

cal  and  meridional  winds  in  the  stratosphere,  (2)  a  would  be 
determined  from  corresponding  wind  data  (w,  v)  in  the  strato¬ 
sphere,  and  (3)  a*  would  be  derivable  from  a  in  the  same  manner 
as  the  eddy  winds.  Because  of  the  unmeasurability  of  the  very 
small  vertical  winds,  they  must  always  be  calculated ,  taking 
into  account  radiation  effects  in  the  stratosphere.  It  must 
thus  be  emphasized  that  any  indirect  method  to  determine  the 
vertical  transports  from  observations  of  the  horizontal  winds 
and  temperature  fields  would  require  the  use  of  assumptions, 
including  those  for  radiation  effects  in  the  stratosphere.  Since 
3-D  GCM  wind  data  already  include  diabatic  effects  in  the  solu¬ 
tion  of  the  primitive  equations,  it  would  be  possible  to  use 
such  wind  data  together  with  Eq.  (14)  so  as  to  derive  estimators 
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of  the  a  and  a*2  statistics  as  a  function  of  latitude,  altitude, 
and  month.  Wind  data  from  the  Geophysical  Fluid  Dynamics  Labora¬ 
tory  (GFDL)/GCM  at  six-hour  intervals  for  each  month  was  available 
for  this  purpose.  The  basic  characteristics  of  this  GCM  have 
been  described  by  Holloway  and  Manabe  (1971),  Manabe  and  Holloway 
(1975),  and  Manabe  and  Mahlman  (1976).  This  GCM  is  characterized 
by  (a)  consideration  of  seasonal  effects;  (b)  a  high  horizontal 
resolution  of  about  265  km,  which  is  Important  for  tracer  trans¬ 
ports  by  wave  numbers  n  >  10  (Mahlman,  1975);  and  (c)  a  vertical 
resolution  of  about  3  km.  from  the  lower  stratosphere  to  the 
middle  troposphere. 

However,  basic  difficulties  in  the  application  of  Eq.  (14) 
were  found  to  be :  (1)  the  fact  that  whenever  v*  decreases 

faster  than  w* ,  the  absolute  magnitude  |a|  of  the  angle  of  the 
mixing  length  path  becomes  relatively  large  and  (2)  the  sta¬ 
tistics  of  a  and  a*2  would  become  then  dominated  by  few  large 
values  of  |a|.  These  results  are  illustrated  in  Tables  5  through 
7  for  low  (24°  N) ,  middle  (48°  N),  and  high  (72°  N)  latitudes, 
respectively.  The  data  in  each  of  these  tables  are  based  on 
sample  GFDL/GCM  wind  data  for  a  6-hour  period  (See  Appendix). 

These  tables  show  the  dominance  of  a  few  large  angles  on  the 
zonal  averages  of  [a]  and  [a*2]  as  a  function  of  pressure  level 
in  the  critical  region  for  the  exchange  of  tracers  between  the 
lower  stratosphere  and  upper  troposphere.  Each  table  shows  (1) 
the  largest  absolute  magnitude  of  a  (or  |a  |)  that  takes  place 
along  discrete  longitudinal  points  for  a  given  latitude  and 
pressure  level;  (2)  the  number  of  longitudinal  points  at  each 
latitude  and  pressure  level,  depending  on  the  magnitude  of  |a|, 
i.e.,  (a)  for  all  magnitudes  including  those  for  |oi  I,  and  (b) 
excluding  angles  that  are  larger  than  either,  say,  5  or  1  deg; 
and  (3)  the  corresponding  [a]  and  [a*2]/[a]2  values  when  includ- 
ing  |  ot  |  or  excluding  (|<x|  <  5,  1  deg)  large  magnitudes  of  |o|. 
Thus,  Table  5  shows  that  at  24°  N  latitude  in  the  tropopause 
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TABLE  5.  SAMPLE  [a],  [a*2]  AT  24°N  LATITUDE 


|«|1iaxl  AND  NUMBER  OF  LONGITUDINAL  POINTS  AS  FUNCTION  OF  |a| 


PRESSURE 

LEVEL,  mbar (km) 

38(22.4) 

65(19) 

110(15.6) 

190(12.1 ) 

|amax  |  , 

deg 

3.1 

14.9 

71  .4 

10.3 

!“!  5  I 

amax  1 

,  Pts 

138 

138 

138 

138 

M  *  5 

deg. 

pts 

138 

133 

135 

137 

|  a  |  £  1 

deg. 

pts 

131 

126 

131 

123 

ZONAL 

AVERAGE  [a](10"4) 

|a|  *  | 

amax  1 

-3.7 

5.0 

99.8 

1  .2 

M  *  5 

deg 

-3.7 

4.5 

7.2 

-11.9 

|  a)  iS  1 

deg 

3.0 

-0.5 

7.9 

-9.2 

_ [a*2 ]/[a32 _ 

| a |  £  |amax|  544.4  5254.0  149.7  26861.1 
|  a |  £  5  deg  544.4  282.5  148.5  1  26.5 
| a |  *  1  deg  250.4  1  1  642.5  35.5  31.4 
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TABLE  6.  SAMPLE  [a],  [ct*2]  AT  48°N  LATITUDE 


let  I  AND  NUMBER  OF  LONGITUDINAL  POINTS  AS  FUNCTION  OF  |a| 
1  max'  _  _ 


PRESSURE  LEVEL 

,  mbar(km) 

38(22.4) 

65(19) 

110(15. 

6)  190(12.1) 

1 amax 1 »  deg 

6.9 

14.9 

6.2 

4.8 

*  1 ama x I ’  pts 

99 

99 

99 

99 

| a |  £  5  deg,  pts 

98 

96 

98 

99 

|a|  £  1  deg,  pts 

89 

94 

96 

93 

ZONAL  AVERAGE  [a](10'4) 


1 « 1 

£ 

1 amax ^ 

-31  . 1 

55.7 

1 

oo 

-2.7 

M 

£ 

5  deg 

-19.2 

-7.7 

2.8 

-2.7 

|a| 

£ 

1  deg 

-0.7 

-8.8 

2.9 

-4.7 

[a* 2 ]/[a  ] 2 


M 

£ 

1 ama x  ^ 

26.2 

46.6 

232.8 

1687.8 

|a| 

£ 

5  deg 

31.4 

63.0 

314.7 

1687.8 

M 

£ 

1  deg 

2477.3 

26.6 

185.4 

34.8 
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TABLE  7.  SAMPLE  [a],  [a*2]  AT  72°N  LATITUDE 


|  ama x  |  AND  NUMBER  OF  LONGITUDINAL  POINTS  AS  FUNCTION  OF  |a| 


PRESSURE  LEVEL,  mbar(km) 

38(22.4) 

65(19) 

110(15.6) 

190(12.1 ) 

1 “max 1 

,  deg 

1.2 

18.1 

0.7 

1  .4 

|a|  * 

1  “max  1 

,  pts 

45 

45 

45 

45 

M  * 

5  deg, 

pts 

45 

44 

45 

45 

M  * 

1  deg, 

pts 

43 

43 

45 

44 

ZONAL  AVERAGE  [a](10'4) 


1  a  1  i 

1  “maxi 

-13.0 

-75.8 

-3.4 

2.2 

1  Ot  I  s 

5  deg 

-13.0 

-7.8 

-3.4 

2.2 

M  * 

1  deg 

-13.9 

-12.6 

-3.4 

-3.4 

[a»2  3/0  ]2 


loti 

s  l“maxl 

19.0 

37.9 

52.4 

360.1 

lot  | 

&  5  deg 

19.0 

23.7 

52.4 

360.1 

i“l 

&  1  deg 

8.0 

2.7 

52.4 

32.8 
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region  (110  mbar  or  15-6  km),  the  largest  absolute  magnitude 
of  a  is  as  high  as  71.4  deg  for  138  discrete  longitudinal  points 
at  this  latitude;  however,  the  bulk  of  the  data  (135  points  or 
98  percent)  had  a  values  smaller  than  5  deg.  The  corresponding 
results  for  [a]  show  that  its  magnitude  drops  from  99.8  x  10"^  to 
7.2  x  10  when  excluding  only  2  percent  of  the  data  with 
angles  |a|  larger  than  5  deg.*  Similar  results  are  obtained 
above  or  below  the  tropopause.  At  65  mbar  (19  km),  the 
[a*2]/[a]2  values  drop  from  5254.0  to  282.5  when  excluding  4 
percent  of  the  data  with  jot)  >5  deg;  whereas  at  190  mbar 
(12.1  km),  [a*2]/[a]2  drops  from  26861.1  to  126.5  when  excluding 
only  1  percent  of  the  data  with  | ot |  >5  deg. 

The  data  in  Tables  6  and  7  show  that  the  number  of  discrete 

longitudinal  points  in  the  numerical  grid  of  the  GCM  for  a  given 

latitude  decreases  as  the  latitude  increases.  Table  6  indicates 

that  (1)  |  ci  I  at  48°  N  and  65  mbar  (19  km)  is  14.9  deg  and 

max  j.  j. 

(2)  the  magnitude  of  [a]  drops  from  55.7  x  10  to  -7.7  x  10“ 
when  excluding  only  3  percent  of  the  data  with  a  >  5  deg.  Table  7 
shows  similar  corresponding  results  at  72°  N  and  the  same  pressure 
level  of  65  mbar  (19  km).  Thus,  the  results  in  Tables  5  through 
7  indicate  that  the  estimators  for  the  a  and  a*2  statistics 
would  be  dominated  by  a  small  percent  of  data  characterized  by 
large  absolute  magnitudes  of  a.  Hence,  the  use  of  the  full  en¬ 
semble  of  120  sets  of  6-hr-period  GCM  wind  data  for  each 
stratospheric  meridional  location  and  month  could  not  be  justified 

The  foregoing  results  in  Tables  5  through  7  are  obtained  from 
a  atriot  use  of  the  Reed  and  German  Eqs .  (6)  for  the  a  statistics 
and  use  of  sample  GCM  wind  data  for  a  6-hour  period,  which 
might  be  somewhat  noisy  due  to  the  rather  low  vertical  resolution. 


Note,  from  Eq.  (14),  that  |a|  depends  on  w  and  v;  which  in  turn 
depends  on  the  number  of  discrete  points  dropped  from  the  data. 

A  computer  code  to  process  the  data  is  also  given  in  the  Appendix 


of  the  GCM  model.  Therefore,  a  determination  of  the  three 
fundamental  parameters  (Kyy,  5,  5*2)  in  the  basic  Reed  and  Ger¬ 
man  formulation  (Eqs.  7  through  10)  would  require  both  reinter¬ 
pretation  of  basic  relevant  concepts  for  the  calculation  of  such 
statistics  and  reformulation  of  the  current  methodology  that 
uses  the  flux  equations  (Eqs.  4  and  5)  for  determining  the 
magnitude  of  the  eddy  transports  coefficients  illustrated  in 
Tables  1  to  3. 
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V.  REED  AND  GERMAN  CONSTRAINTS  CARRIED  OVER  TO  TWO-DIMENSIONAL 

PHOTOCHEMICAL  MODELS 


A  need  for  reformulation  of  the  current  methodology  for 
establishing  both  the  magnitude  and  even  nature  of  2-D  eddy 
transport  coefficients  is  emphasized  by  the  following: 

(1)  Even  if  the  a  and  oi* 2  statistics  were  determin¬ 
able  from  wind  data  through  Eq .  (6),  the  system  of 

Eqs.  (7)  through  (10)  would  then  specify  K  _  and  K 

y  z  zz 

in  terms  of  K  ;  which  could  be  determined  from 

«/  i/ 

Eq.  (4).  The  result  would  be  an  overspecified  system 
of  two  equations  (A  and  5)  for  one  unknown  (K  )  in- 

J  tf 

stead  of  the  undetermined  system  of  two  equations 

(4  and  5)  and  three  unknowns  (K„,r,  K.  K„_)  ,  as 

yy  yz 

dealt  with  in  the  current  methodology  for  the  2-D 
parameterizatlons  of  tracer  transports  in  Tables  1 
through  3. 

(2)  As  shown  by  Andrews  and  McIntyre  (1976), 

Eq.  (3)  is  valid  for  waves  with  small  amplitudes. 

However,  as  pointed  out  by  Holton  (1980),  Z and  Zz 
are  to  be  interpreted  as  fluid  particle  displacements 
caused  by  the  waves.  Hence,  if  v*  and  w*  denote  the 
wave  velocity  components 


.ffz 

dt 


w*  = 


d  Z 

z 

dt 


(15) 
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Then,  in  place  of  Eqs.  (7)  through  (10)  the  follow¬ 
ing  result  is  obtained 


K  =  v*£  =  L  —2.  =  —  ■  y 

yy  y  y  dt  2  dt 


(16) 


d  1 

K  =  v*T"  -  l  J 
yz  z  z  dt 


(17) 


Kzy  -  -  l 


K  _  =  w»T“  =  l 


-ax 

y  St¬ 

ar 


z  d 


(£  £  ) 
dt  y  z 


.  d£* 
1  z 


-  K 


yz 


zz  z  z  dt  2  dt 


(18) 


(19) 


Equation  (18)  shows  that,  in  general,  K  ^  K  ;  a 

zy  y  z 

result  that  is  different  from  the  Reed  and  German 
Eq.  (9),  which  is  derived  from  the  critical  assumption 
that  a*  is  independent  of  v  and  £. 

The  foregoing  results  indicate  that  the  reform¬ 
ulation  of  the  methodology  to  derive  2-D  eddy  co¬ 
efficients  may  further  involve  the  determination  of 
the  four  parameters  given  by  Eqs.  (16)  through  (19), 
instead  of  the  three  in  the  Reed  and  German  Eqs.  (7) 
through  (10).  These  results  indicate  again  that  the 
current  2-D  diagnostic  parameterizations  of  strato¬ 
spheric  tracer  transports  (as  illustrated  in  Tables 
l  through  Z),  based  on  Eqs.  (7)  through  (10)  or  Eq.  (9) 
with  (4),  (5),  and  (2)  for  either  inert  or  chemically 
active  tracers ,  may  not  rely  on  the  parent  Reed  and  Ger¬ 
man  (1965)  formulation  to  justify  their  prognostic 
use  in  the  forecasting  of  adverse  anthropogenic 
effects  on  atmospheric  ozone. 
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VI.  PAST  INTERPRETATIONS  OF  REED-GERMAN  RESULTS 


Previous  attempts  (Mahlman,  1975;  Hidalgo,  1978)  to  assess 
the  prognostic:  validity  of  diffusivity  concepts  used  in  2-D 
parameterizations  of  tracer  transports  in  the  stratosphere  have 
been  based  on  GFDL  GCM/tracer  model  data  for  a  tracer  character¬ 
ized  as  being  chemically  inert  in  the  stratosphere  but  with 
heterogeneous  sinks  in  the  troposphere.  The  GCM/tracer  model 
combines  the  GCM  numerical  solutions  from  the  primitive  equa¬ 
tions  with  those  from  the  3-D  continuity  Eq.  (1)  for  a 
chemically  inert  tracer,  i.e.,  the  auxiliary  tracer  model 
solves  for  the  time  evolutions  of  tracers  from  dR/dt  =  0  while 
using  as  inputs  the  wind  field  of  the  parent  GCM.  Because  of 
constraints  on  computer  resources,  the  same  annual  cycle  from 
the  GCM  is  reused  in  subsequent  cycles  of  the  tracer  model. 

While  this  procedure  does  not  allow  use  of  the  annual  variability 
of  the  wind  fields,  the  annual  variability  of  the  tracer  is  soley 
due  to  the  tracer  itself.  The  GCM/tracer  model  data  used  to  assess 
the  prognostic  validity  of  diffusivity  concepts  in  2-D  para¬ 
meterizations  have  then  come  from  the  following  three  numerical 
experiments : 

(1)  A  "vertical  stratification"  experiment,  which  utilized 
a  tracer-source  in  the  upper  stratosphere.  This  was  the 
earliest  ozone-precursor  numerical  experiment  that  was  de¬ 
signed  to  provide  a  crude  simulation  of  the  equilibrium 
photochemical  stratification  of  ozone.  It  specified  a 
constant  mixing  ratio  in  the  upper  stratosphere  of  a 
chemically  inert  tracer,  which  was  assumed  to  have  hetero¬ 
geneous  sinks  in  the  lower  troposphere  for  its  removal  by 
rainout  processes  (Mahlman,  1973a,  1975). 
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(2)  A  "raid-latitude  point  source"  experiment,  which 
utilized  an  instantaneous  release  of  a  chemically  inert 
tracer  at  a  point  located  at  a  middle  northern  latitude 
in  the  lower  stratosphere.  This  experiment  simulated  the 
distributions  of  radioactive  nuclear  debris  in  the  strato¬ 
sphere  from  a  nuclear  explosion.  It  specified  an 
initial  distribution  of  mixing  ratio  at  and  near  the  point- 
source  of  a  tracer,  which  was  again  characterized  as 
being  chemically  inert  in  the  stratosphere  but  with  hetero¬ 
geneous  sinks  in  the  troposphere.  (Mahlman,  1973a,  1975; 
Mahlman  and  Moxim,  1978). 

(3)  An  SST  (supersonic  transport )-like  experiment,  which 
utilized  constant-line  sources  of  a  chemically  inert  tracer 
located  at  specified  zonal  segments  at  a  northern  middle 
latitude  in  the  lower  stratosphere  (Mahlman,  1973b,  1978; 
Hidalgo,  1978).  This  experiment  simulated  the  time  (years) 
evolutions  of  the  monthly  distribution  of,  say,  N0y 

(NO  +  N02  +  2N20^  +  NO^  +  HNO^)  as  the  result  of  emis¬ 
sions  of  SST  engine  effluents  in  the  lower  stratosphere. 
Again,  the  characteristics  of  the  inert  tracers  are  as 
in  the  foregoing  two  experiments. 

The  approach  used  in  these  assessments  of  the  prognostic 
validity  of  2-D  parameterizations  of  eddy  transports  consisted 
in  the  use  of  the  Reed  and  German  Eqs .  (4),  (5),  and  (9),  to¬ 
gether  with  3-D  GCM/tracer  model  data  from  the  above  numerical 
experiments.  The  3-D  tracer  data  was  averaged  using  the  (~) 
operator  to  obtain  estimators  for  the  2-D  tracer  statistics  of 
the  meridional  and  vertical  eddy  transports  and  corresponding 
gradients  of  mixing  ratio.  However,  this  use  of  the  Reed  and 
German  equations  with  averaged  3-D  GCM/tracer  model  data  could 
not  overcome  the  difficulty  of  defining  properly  the  magnitude 
of  the  ci  and  a*2  statistics;  i.e.,  the  problems  faced  previously 
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when  attempting  to  solve  the  undetermined  system  of  the  two 
eddy  flux  equations  for  the  three  unknown  eddy  coefficients. 


Since  Eq.  (14)  yields  a  =  w#/v#  for  small  angles  of  a,  an 
approximation  used  to  estimate  a.  was  to  multiply  the  numerator 
and  denominator  of  a  by  R*  and  assume  that  a  =  w#R#/v*R*'.  Use 
of  Eqs.  (4),  (5),  and  (9)  yielded  then  (Mailman,  1975): 


v^R*  p  +  w*R*  ~ 

9y  9z 

and 


K 


zz 


W*R#  2 


^FrF 


3R  .  TJWuW  3R 

__  +  w*R* 


(21) 


The  use  of  Eqs.  (20),  (21)  and  averaged  3-D  GCM/tracer  lata  for 
the  eddy  transports  and  gradients  of  mixing  ratio  in  t-hese 
equations  would  then  yield  the  magnitude  of  the  eddy  transport 
coefficients.  Such  results,  based  on  the  "vertical  stratifica¬ 
tion"  and  "mid- latitude  point  source"  experiments,  indicated 
that  there  were  regions  in  the  stratosphere  where  the  eddy 
transport  coefficients  became  even  negative  (Mahlman,  1975); 
results  that  were  also  confirmed  for  both  the  mid-latitude  or 
instantaneous-point  source  as  well  as  the  constant-line  or  SST 
sources  as  indicated  in  Tables  8  through  10  (Hidalgo,  1978). 
However,  these  latter  results  also  indicated  that  the  correspond¬ 
ing  magnitudes  of  a  became  different  for  the  point  source  and 
constant-line  source  experiments  (Table  8)  at  the  same  strato¬ 
spheric  location  and  month  (season)  instead  of  remaining  fixed 
by  the  variability  of  the  stratospheric  winds  as  by  Eq .  (14). 
Furthermore,  since  division  of  Eq .  (21)  by  (20)  yields 
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8.  SLOPE  OP  «»??#L|*TH(JJI!j)JS  8,VEN  8V  PP.MUm 
Source:  Mahlman,  Private  Communication .  I97fi 
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Por  instantaneous  Injection  at  38  mbar  and  72°N,  for  example,  [«] 


TABLE  9.  VERTICAL  EDDY  COEFFICIENTS  AS  GIVEN  BY  PRIMITIVE 
EQUATIONS,  Kzz  ( 10*cm2/sec)* 

(Source:  Mahlman,  Private  Communication,  1978) 
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1. 


K  /K  =  a2,  the  use  of  Eqs.  (7)  and  (8)  gives  then  K  K  /K2  = 

zz  yy  °  y.y  zz  yz 

Hence,  Eq.  (13)  indicates  that  the  assumption  a  =  w*W/v*W  leads 
to  a*2  =  0.  Since  the  use  of  Eqs.  (20)  and  (21)  contradicts 
both  the  sole  dependency  of  a  on  the  wind  variability  as  well 
as  the  apparent  undeterminability  of  a*2  from  wind  data,  pre¬ 
vious  conclusions  of  negative  eddy  transport  coefficients  based 
on  these  equations  are  not  warranted  from  a  strict  use  of  the 
Reed  and  German  formulation. 

It  should  however  be  noted  that  the  foregoing  implied 
assumption  of  a*2  =  0  in  Eqs.  (20)  and  (21),  as  well  as  the  re¬ 
sults  of  Tables  8  through  10  is  as  in  the  Crutzen  model  at  middle 
and  high  latitudes  (Table  1).  On  this  basis,  it  appears  that 
the  use  of  a*2  =  0  in  diagnostic  parameterizations  (for  a  strato¬ 
sphere  without  SST  pollutants)  would  likewise  yield  a  values 
that  are  inconsistent  with  those  required  by  the  Reed  and  German 
equations  for  a  stratosphere  with  SST  pollutants  (Table  8); 
i.e.,  such  diagnostic  parameterizations  with  fixed  a  values  at 
a  given  stratospheric  location  and  month  (season)  may  have  a 
questionable  prognostic  value.  Furthermore,  the  validity  of 
even  the  diagnostic  value  of  2-D  parameterizations  based  cn 
a*2  =  0  are  open  to  question,  since  the  foregoing  use  of  GCM/ 
tracer  model  data  together  with  the  Reed  and  German  equations 
requires  that  the  eddy  coefficients  may  become  negative  (Tables  9, 
10).  Moreover,  it  is  of  interest  to  Indicate  that  GCM/tracer 
model  data  suggest  that  the  magnitude  of  the  tracer  eddy  trans¬ 
ports  may  depend  on  the  prevailing  time  and/or  space  tracer 
configuration  in  the  stratosphere  (Mahlman  and  Moxim,  1978; 
Hidalgo,  1978);  a  result  that  would  affect  the  accuracy  of  diag¬ 
nostic  2-D  parameterizations  (Tables  2,  3)  that  are  not  based 
on  a*2  *  0. 
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VII.  OUTLOOK  FOR  PROGRESS  IN  PROGNOSTIC  2-D  FORMULATIONS  OF 

TRACER  TRANSPORTS 

Recent  developments  stemming  from  a  Lagrangian  instead  of 
a  Eulerian  description  of  wave  motions  and  associated  mean 
flows  (Andrews  and  McIntyre,  1978  a  and  b;  Dunkerton,  1978)  indicate 
that  (a)  the  net  countergradient  "diffusion"  is  really  an  ad- 
vection  by  the  Stoke's  drift  due  to  the  waves  (Wallace,  1978); 

(b)  the  entire  bulk  motion  of  the  center  of  mass  of  air  parcels 
and  tracers  would  be  described  advectively,  not  by  Kyz  slantwise 
diffusion;  however  there  is  still  need  for  diffusion  to  explain 
the  lateral  and  vertical  spread  of  a  tracer  about  its  center  of 
mass  (Holton,  1980);  and  (c)  the  Reed  and  German  form  of  Eqs. 

(4)  and  (5)  might  be  retainable  for  small-amplitude  waves. 

However,  a  rather  drastic  reinterpretation  of  the  physical  mean¬ 
ing  of  Eqs.  (4)  and  (5)  becomes  necessary  (Plumb,  1979;  Matsuno, 
1980);  namely,  major  parts  of  the  meridional  eddy  fluxes  of 
ozone,  potential  temperature,  and  other  vertically  stratified 
quantities  would  be  the  result  of  a  systematic  correlation  be¬ 
tween  the  vertical  displacement  and  the  meridional  velocity  of 
air  parcels  arising  from  a  particular  structure  of  planetary 
waves,  and  hence  the  fluxes  are  nearly  transverse-gradient 
rather  than  obliquely  down-gradient  (Clark  and  Rogers,  1978). 

Such  eddy  transports  are  shown  to  be  nearly  advective  instead 
of  diffusive  in  nature  and  may  be  expressed  by  eddy  transport 
coefficients  with  Kyz  /  Kzy- 

The  foregoing  possibility  is  suggested  from  exploratory 
investigations  that  considered  the  stratospheric  motions  in¬ 
duced  by  an  upward,  steady-state  propagation  of  a  planetary 
wave  from  the  troposphere  into  a  stratosphere  characterized  by 
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a  uniform  zonal  (i.e.  u  or  eastward)  circulation  in  a  latitude 
domain  (Matsuno,  1980).  The  approach  consisted  in  obtaining 
(1)  solutions  of  the  potential  vorticity  equation  for  the  wave 
structure  and  the  induced  motions  in  the  stratosphere,  (2)  La- 
grangian  trajectories  of  the  air  parcel  motions  and  (3)  eddy 
transports  of  x  (e.g.  ozone)  carried  by  the  air  parcels  in  a 
form  which  is  proportional  to  local  gradients  of  x  as  in  Eqs. 

(4)  and  (5).  The  important  relevant  results  from  this  mechan¬ 
istic  approach  are  then  as  follows:  (1)  the  Lagrangian  tra¬ 
jectories  of  air  parcels  consists  of  cyclic  motions  along 
ellipses  as  well  as  the  linear  motions  (as  shown  in  Pig.  1)  in 
the  meridional  domain  of  the  stratosphere,  (2)  the  orientation 
of  the  major  axis  of  the  meridional  ellipses  is  indeed  as  that 
of  a  in  Pig.  1  and  with  the  same  order  of  magnitude  (~10~4)  as 
in  the  Reed  and  German  concept,  (3)  if  the  meridional  ellipses 
were  superimposed  on  Pig.  1,  the  direction  of  motion  of  air 
parcels  would  be  clockwise  at  any  latitude,  (4)  the  eddy  trans¬ 
port  coefficients  may  be  expressed  in  terms  of  gradients  as  in 

Eqs.  (4)  and  (5)  but  with  unequal  K  and  K  coefficients  and 

y  “  z  y 

(5)  the  eddy  transport  coefficients  are  a  function  of  the  period 
of  the  cyclic  motions  along  the  ellipses  and  the  time  constant 
of  chemical  adjustment.  Hence,  these  encouraging  exploratory 
results  suggest  that  a  prognostic  2-D  formulation  of  eddy  trans¬ 
ports  of  chemically  inert  or  active  tracers  in  the  stratosphere 
might  be  feasible. 
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•VIII.  TESTS  OF  1-0  PARAMETERIZATION 


The  1-D  photochemical  models  of  the  stratosphere  and  tropo¬ 
sphere  have  played  a  dominant  role  during  recent  years  in  crude 
assessments  of  ozone  effects  from  emissions  of  (1)  aircraft  engine 
effluents  at  high  altitudes  [Grobecker  et  al.,  197^;  National 
Academy  of  Sciences  (NAS),  1975;  Oliver  et  al.,  1977;  Broderick, 
1977]  and  (2)  chlorofluoromethanes  (CFMs)  at  the  earth's  surface 
(NAS,  1976;  NASA,  1977;  Crutzen  et  al . ,  1977).  Basic  factors  for 
the  popularity  of  1-D  photochemical  models  have  been  the  uncertain 
knowledge  of  the  relevant  chemistry  together  with  their  capability 
for  a  wide  coverage  of  chemical  cycles  and  species  as  a  result  of 
their  rather  modest  requirements  for  computer  and  manpower  re¬ 
sources  . 

The  1-D  continuity  equations  for  the  conservation  of 
chemical  species  are  obtained  by  taking  the  global-annual  sta¬ 
tistical  average  of  Eqs.  (1)  for  an  ensemble  of  many  years. 

If  { (  ) }  denotes  the  global-annual  ensemble  average  of  (  ) , 
then  {(  )}  represents  a  generalization  of  the  2-D  ensemble 
average  (“)  by  extending  it  with  respect  to  latitude  and  even 
time  from  a  monthly  to  an  annual  basis.  Although  it  would  be 
possible  to  likewise  define  (  )**  =()-{()}  so  as  to  separ¬ 
ate  at  each  altitude  the  1-D  average  of  the  vertical  eddy  trans¬ 
ports  of  mixing  ratio  {w**R**}  from  those  of  the  corresponding 
vertical  winds  or  1-D  circulation  (w}{R),  it  is  customary  to 
extend  the  concept  of  2-D  eddy  transports  to  the  1-D  case  by 
assuming  that  the  1-D  transport  coefficient  Kz  (or  K  )  includes 
now  the  transports  by  the  1-D  vertical  circulation.  Since 
{^p*R^U}  =  0,  where  5  denotes  the  local  horizontal  wind  vector, 
the  basic  continuity  equations  for  1-D  photochemical  models  in 
pressure  coordinates  is  given  by 

39 


at  ^Ri*  +  ap  *wRi*  =  ~ 

where  p  denotes  pressure  and  w  =  dp/dt .  Again,  with 

{uRi}  =  -  Kp  |p  {Ri}  *  (23) 

where  Kp  is  given  in  (mbar)Vsec,  Eq.  (22)  becomes  a  set  of 
differential  equations  for  the  mixing  ratio  Rj.  The  above 
equations  can  be  written  in  the  traditional  altitude  instead  of 
pressure  coordinate  through  use  of  the  hydrostatic  equation. 

The  1-D  result  equivalent  to  that  of  Eq.  (2)  for  a  chemically 
inert  tracer  in  the  stratosphere  takes  then  the  following  form: 

It  {R}  =  f  fz  5  {wR}  »  (24) 

where  p  is  proportional  to  {p/T}.  Equation  (23)  becomes 

{wR}  -  -  Kz  {R}  .  (25) 

Since  1-D  photochemical  models  assume  that  the  1-D  transport  of 
chemical  species  is  given  by  Eq.  (25)  for  a  chemically  inert 
tracer,  it  is  seen  that  the  usual  definition  of  the  1-D  trans¬ 
port  coefficient  K_(z)  is  given  as  a  ratio  of  the  global-annual 
ensemble  average  of  the  flux  of  mixing  ratio  by  both  the  vertical 
components  of  the  circulation  and  eddy  motions  to  the  vertical 
gradient  of  the  corresponding  average  of  mixing  ratio.  The 
lack  of  worldwide  observations  for  an  ensemble  of  many  annual 
cycles  of  a  suitable  tracer  in  the  critical  regions  of  the  lower 
stratosphere  and  upper  troposphere  has  made  it  very  difficult 
to  arrive  at  a  generally  acceptable  K  (z)  profile  for  the  fore- 
casting  of  long-range  adverse  anthropogenic  1-D  ozone  effects. 


HO 
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The  continuity  equations  of  1-D  photochemical  models  have 
the  form  of  Eqs.  (24)  and  (25),  except  that  replaces  R  and 
the  right-hand  side  of  Eq.  (24)  must  incorporate  the  production 
and  loss  terms  {Pi>  and  {L^}.  A  proper  physical  interpretation 
of  1-D  ozone  effects  is  then  very  difficult,  because  such  re¬ 
sults  are  based  on  global-annual  ensemble  averages  of  not  only 
the  atmospheric  transports  but  also  of  the  photochemistry  and 
shortwave  solar  radiation  (or  solar  zenith  angle),  chemical 
kinetics,  heterogeneous  chemistry  and  emissions  of  aircraft 
engine  effluents.  However,  the  common  practice  of  using  1-D 
photochemical  models  to  forecast  anthropogenic  ozone  effects  mo¬ 
tivates  interest  in  exploring  the  characteristics  of  1-D  transport 
coefficients  in  the  critical  region  of  the  lower  stratosphere  and 
upper  troposphere  for  such  effects.  This  can  again  be  done  by 
using  3-D  GCM/tracer  model  data,  i.e.,  by  obtaining  estimators 
from  the  1-D  global-annual  average  of  3-D  GCM/tracer  data  for 
the  vertical  fluxes  and  vertical  gradients  of  mixing  ratio  as  re¬ 
quired  by,  say,  Eq.  (23)  with  R.^  =  R  for  a  chemically  inert 
tracer. 

Recent  results  from  numerical  experiments  at  the  Geophysical 
Fluid  Dynamics  Laboratory  allow  a  rather  comprehensive  display 
of  the  characteristics  of  1-D  transport  coefficients  as  a  func¬ 
tion  of  (1)  tracer  configuration  in  the  stratosphere  with  regard 
to  both  the  space  and  time  (years)  domains. and  (2)  tracer-source 
location  in  the  stratosphere  as  well  as  at  the  earth’s  surface. 

In  addition  to  the  three  numerical  experiments  described  pre¬ 
viously,  there  are  two  more  recent  experiments.  One  is  called 
"simple  ozone",  which  is  a  refinement  of  the  "vertical  stratifi¬ 
cation"  experiment  for  a  tracer-source  in  the  upper  stratosphere. 
In  this  subsequent  experiment,  use  was  made  of  "simple"  chemistry 
at  the  top  of  the  model  together  with  an  instantaneous  relaxation 
to  a  jpecified  observed  average  ozone  value  at  10  mbar,  but  with 
otherwise  the  same  inert  and  heterogeneous  tracer  characteristics 
as  in  the  "vertical  stratification"  experiment  (Mahlman,  Levy, 
and  Moxim,  1979).  The  other  recent  numerical  experiment  used  a 


small  uniform  source  of  NgO  (15  M  ton/year)  at  the  earth's 
surface  (Levy,  Mahlman,  and  Moxim,  1979)  and  is  of  direct  inter¬ 
est  to  the  CFM  problem;  because  as  CFMs,  N20  propagates  upwards 
from  the  earth's  surface  to  the  stratosphere  after  a  long  resi¬ 
dence  time  (several  decades)  in  the  troposphere. 

Some  characteristics  of  monthly-annual  estimators  for  1-D 
transport  coefficients  based  on  early  "stratification"  and 
"mid- latitude  point  source"  results  have  been  given  by  Mahlman 
(1975).  These  results  Indicated  that  the  K  (z)  profiles  were 
not  unique,  i.e.,  the  1-D  transport  coefficient  was  found  to  be 
strongly  dependent  upon  the  tracer  configuration  in  the  strato¬ 
sphere.  Furthermore,  these  results  indicated  that  the  Kz(z) 
profile  may  not  necessarily  be  independent  of  time  (years)  as 
assumed  in  1-D  photochemical  models  during  th«.  transition  years 
between  the  equilibrium  states  of  an  initial  'polluted  to 
polluted  stratosphere. 

Table  11  shows  the  1-D  transport  coefficients  derived  from 
the  GFDL  3-D  GCM/tracer  model  experiments  as  a  function  of 


TABLE  11.  1-D  TRANSPORT  COEFFICIENTS  Kz(z),  lOV/sec  AS  A  FUNCTION  OF 

SOURCE  LOCATION  AND  TIME  (YEARS)  IN  THE  LOWER  STRATOSPHERE 
AND  TROPOSPHERE  (Mahlman,  1980) 


LEVEL  ALTITUDE* 

POINT-SOURCE 

LINE- 

SOURCE 

STRATIFICATION 

"SIMPLE  03" 

n2o 

(mbar) 

(km) 

Year  1 

Year  4 

Year  7 

Year  9 

Year  1 

Year  4 

Year  4 

27.6 

24.2 

8.44 

1.49 

2.43 

-2.16 

0.52 

0.74 

0.60 

1.21 

52.3 

20.2 

0.10 

0.77 

1.99 

-2.36 

0.42 

0.61 

0.66 

0.62 

80.7 

17.6 

0.65 

1.26 

1.26 

1.41 

0.92 

1.10 

1.14 

1.01 

149.9 

13.7 

1.92 

1.25 

1.86 

1.92 

1.52 

1.11 

1.16 

0.84 

240.6 

10.5 

7.02 

3.67 

12.94 

13.00 

5.19 

3.91 

3.94 

10.77 

412.0 

6.8 

9.53 

6.24 

11.54 

11.61 

12.53 

6.00 

6.13 

9.57 

606.7 

4.0 

10.04 

11.29 

44.03 

42.34 

18.31 

12.32 

12.80 

_  ? 

Standard  altitude. 

★* 

Equilibrium  conditions  (reached  after  -  2000  years). 


tracer-source  location  and  tracer  configuration,  i.e.,  tracer- 
source  characteristics  and/or  time  (years).  The  results  in 
'ble  11  indicate  the  following: 

•  For  the  point  and  line  sources  in  the  lower  stratosphere 
of  interest  for  aircraft  ozone  effects,  the  1-D  trans¬ 
port  coefficients  for  the  lower  stratosphere  can  be  a 
function  of  both  source  characteristics  and  year.  For 
example,  at  20.2  km,  the  two  instantaneous-point  source 
columns  show  that  the  value  of  the  1-D  transport  co¬ 
efficient  for  the  first  year  can  increase  by  a  factor 

of  about  eight  in  the  subsequent  year;  whereas  the  corres¬ 
ponding  constant-line  source  columns  indicate  that  the 
1-D  transport  coefficient  for  the  ninth  year  can  become 
negative,  a  fact  that  implies  distortions  of  the  K(z) 
profile  (Mahlman,  1980). 

•  The  altitude  for  the  minimum  value  of  the  1-D  transport 
coefficient,  which  is  of  critical  importance  for  pre¬ 
dictions  of  anthropogenic  ozone  effects,  can  be  a  function 
of  source  location.  For  example,  the  N20  column  indicates 

that  the  altitude  for  minimum  K  for  N~0  is  at  about 

z  2 

13.7  km;  which  is  lower  than  that  (20.2  km)  for  the 
instantaneous-point  source  at  mid-latitudes  in  the  lower 
stratosphere . 

Although  the  results  in  Table  11  are  subject  to  constraints 
that  may  be  introduced  by  the  rather  low  vertical  resolution 
(3  km)  of  the  GCM/tracer  model  in  the  critical  tropopause  and 
lower  stratosphere  regions,  the  trends  established  by  such  ex¬ 
tensive  calculations  may  not  be  ignored  in  the  absence  of  global 
observations  involving  suitable  tracers  for  an  ensemble  of  many 
years.  Hence,  1-D  transport  coefficients  that  are  assumed  to 
be  independent  of  tracer  configuration  (i.e.,  tracer-source  loca¬ 
tion  and  characteristics  as  well  as  years)  as  usually  done  in 
1-D  photochemical  models,  cannot  avoid  serious  distortions  in 
their  prognostic  descriptions  of  transports. 
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APPENDIX 


This  Appendix  provides  the  sample  GFDL/GCM  data  for  a 
6-hour  period  as  a  function  of  latitude,  pressure  level,  and 
longitude.  The  data  shows  the  vertical  (w),  northward  (v), 
and  eastward  (u)  winds  in  cm/sec  at  28°N,  48°N,  and  72°N 
latitude.  Note  that  the  number  of  longitudinal  points  decrease 
from  138  at  the  lowest  latitude  to  99  and  45  at  the  middle  and 
high  latitudes,  respectively.  This  Appendix  also  gives  a  com¬ 
puter  program  for  the  calculation  of  the  [a]  and  [a*2]  statistics 
based  on  both  the  Reed  and  German  definition  (Eq.  14)  and  by 
replacing  arbitrarily  the  northward  wind  by  the  eastward  wind. 

The  former  slope  is  denoted  by  a^,  and  the  latter  by  a^.  The 
purpose  of  computing  a2  was  to  investigate  if  the  [a^]  and 
[a*2]  undeterminability  of  the  statistics  would  be  modified 
from  possible  consideration  of  the  zonal  circulation.  The 
results  were  negative,  i.e.,  the  undeterminability  of  [o^]  and 
[cu*]  was  as  that  of  the  Reed  and  German  [cx-^]  and  [a^*2]. 
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TABLE  A- 1 .  W-VERTICAL  WIND,  LATITUDE  24  DEGREES  NORTH 

PRESSURE  LEVEL 


30MR 

65MB 

11 0MB 

190MB 

-1,1^2 

■3  967 

,7o0 

1.270 

-1 ,077 

:iB« 

>59 

1.149 

.219 

.286 

-.3*0 

.814 

,099 

-,034 

-.919 

1.470 

,823 

.1.230 

-2.198 

.686 

.091 

-1.0*1 

-2.430 

,515 

4.144 

.>99 

-5.157 

1.303 

.1.201 

-.475 

-1.131 

,999 

-4.17? 

-»7ol 

1.663 

,960 

-1.141 

-1.976 

.602 

,475 

,641 

.1.039 

-.965 

1.297 

1.926 

.1,600 

-3.593 

,573 

1,970 

.,206 

-.12® 

-.570 

1.91® 

-.047 

-.202 

-1,634 

-.799 

1.341 

,697 

.1,030 

.1.23? 

2  >38 

1.300 

,2«0 

-1.043 

.1,2*6 

-.108 

2.005 

-1.098 

-1.7U 

.674 

2,679 

.95? 

-2.288 

-.146 

1.470 

-.904 

.090 

-2.017 

.277 

-.001 

1.974 

-1.325 

1.197 

-1.9 Bt 

-.0*6 

3.635 

.777 

-1.674 

-.021 

5.422 

-.120 

-,l7t 

,761 

3.304 

.006 

8.490 

-1.0U7 

-1.397 

•1.072 

4.33* 

>00 

-2.327 

.1.050 

1.190 

3,014 

2.796 

-.274 

•1.043 

1,069 

4.367 

-.096 

1.424 

.lit 

3.267 

.415 

-.soi 

1.919 

2.997 

.723 

.914 

-,  1 03 

-.720 

-.360 

4,69? 

•2,106 

-3.094 

-2.395 

.34? 

4,200 

-1.397 

-.970 

-4.27« 

3,240 

3.741 

,029 

-2.76? 

*1  ,*42 

.003 

.619 

-4.034 

-3, 101 

-3.314 

-.680 

.1,11® 

-1,931 

.1.010 

.1.077 

,634 

2,000 

3,105 

-.665 

-1.40? 

-,7»3 

-1.501 

-,609 

-1.99? 

-.014 

,692 

.,043 

.1.13? 

1 .9*9 

5.700 

.590 

-1.470 

1,202 

2.972 

-.129 

-2.198 

1,260 

4.057 

.019 

-.021 

.2.637 

7.059 

-.214 

•  f  46  J 

1.933 

4.670 

•2,063 

-39404 

4.060 

2.910 
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TABLE  A- 1 .  CONTINUED 


38*S 

-.85) 

.657 

-.240 

.823 

.047 

-.590 

1.340 

.844 

.738 

.579 

.472 

.5*7 

1.0)5 

1.606 

.799 

1.413 

.666 

.218 

-.189 

-1.571 

-2.6?8 

-2.489 

-1.793 

-.906 

.524 

1.295 

.597 

-.403 

-.375 

-.7)7 

-.678 

.398 

.825 

-1.743 

.483 

-1.911 

•1.886 

-.045 

-.4.79 

-.716 

-1.587 

-1.287 

.229 

-.154 

.617 

.260 

-.679 


9588 
-•  916 
1  .447 
-1.052 
-2.639 
-1.00? 
.39* 

1.151 

-.866 

“4.0*? 

-4.124 

.034 

•1.640 

.33! 

1.652 

*.554 

-.024 

.22* 

1.706 

2.701 

.67? 

-4,713 
-4.412 
-4.773 
-4 ,938 
•3.06u 
-4 .846 
1.73S 
>3.947 

-1.552 

.060 

-.388 

-,66i 

.974 

-4,9l« 

-.533 

.005 

.264 
-4.579 
,874 
,821 
-1 ,630 
-.32? 
•1.221 
-.932 
4.344 


11096 
.2/3 
-.933 
2,699 
2.404 
2,127 
2.138 
.6*4 
2.098 
1.605 
2.355 
•  »6o3 
.030 
4,122 
3,828 
1,626 
1  ,«38 
2.096 
1,692 
1,766 
2,168 
*398 
-.241 
-2,228 
-2,196 
-2,4il 
•3,8*7 
-2,5vz 
-2,899 
-4,6*8 
-2,863 
-1,210 
.,237 
•3,469 
.563 
»,6o3 
.064 

2.3*1 

-»0o2 

,«7 

-.840 

,239 

,140 

-1,364 

-.540 

-.004 

-.538 

,168 


190Hti 
—  •  676 
-5.222 
-1.624 
-.428 
-.24) 
-2.846 
-2.114 
.242 
4.270 
6,069 
2.798 
•5,305 
-2.499 
-1.521 
..849 
-1.548 
-1.136 
-2.229 
-1.846 
-1.098 
3.444 
6.545 
2.387 
-1.328 
-5.179 
-3.213 
•  .444 
-4 .266 
-6.360 
-3.737 
-2,634 
-1,783 
-3.310 
-.116 
-4.364 
.992 
.120 
-4.733 
-.332 
-.009 
.,901 
-1.993 
2.125 
1.331 
-1 .844 
-1.314 
•1 .268 
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TABLE  A-l .  CONTINUED 


304* 

•  000 
•1.4?? 
••2*4 
-.043 
•  200 
1.340 
.495 
••300 
••040 
••40? 

.079 
•  .2*0 
.540 

•  .303 

•  .4*9 
1.100 

.592 

.044 

.03? 

«.2(f 

1.033 

-.315 

•1.510 

-.441 

-.001 

1.793 

2.2?4 

.3*1 

-.41? 

,577 

-.173 

1.299 

•1.077 

-.074 

•374 

••4*0 

•575 

.1*7 

«3l* 

•00? 

•074 

••463 

•4*0 

-.112 

•*44? 


65Mb 

•  •240 

•1.189 

•  •080 
2.0*6 
1.851 

.692 

1.245 

2.310 

.109 

•1.39* 

..53S 

.865 

.**« 

.,675 

-.604 

2.325 

•a.»i“ 

..62] 

•1.194 

1.2?9 

3.076 

3.059 

*•5?? 

•4,8i? 

•2.2?] 

3.0*4 

.1.0*1 

•  ,760 
•2.531 
•1.9*4 

a, 71k 

•  ,9|0 
•1.161 

.073 

■••000 

1*519 

.530 

.96? 

«.32* 

.14? 

2,643 

■1.019 

•a. 004 


110*0 
•  ,012 
.595 
1.020 
1.3*2 
.if  4 

.230 

■*2o3 

•1.6*3 

.10* 

1.201 

.062 

.7*2 

-,2s9 

.116 

-,2f4 

-.443 

•  .043 
-.710 
1.019 
1.740 

•1 ,402 

-.710 

.**0 

1,1*9 

-2.043 

-♦.*?2 

—  ,  020 
200 
-.700 
•2.6*1 
.134 
-.303 
-,9t* 
•  *27 
-.2*4 
•.1*7 
•.3  J4 

•  •*45 
•3r* 

1.062 

.390 

•1.273 


190*0 
-.474 
-.946 
••02* 
•9.059 
3?505 
•  .  066 
•5.60? 

.115 

-2.242 
.386 
1.39* 
1.055 
2i456 
7.053 
3.3*5 
1.101 
7.756 
-7.17? 
•7,000 
•1 ,001 
7.994 
7.901 
•1,002 
•6.7*3 
•3.9*2 
6.33? 
1.0*0 
.3.400 
-.291 
-.033 
2.62* 
1.043 
.614 
-.750 
•  26? 

•  *165 
-.643 

•1.255 

•  •194 
•2*110 

•  .13? 
-.77* 

•1.93? 
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TABLE  A- 2 •  V-NORTHWARD  WIND,  LATITUDE  24  DEGREES  NORTH 


PRESSURE  LEVEL 


>329,p46 
>401  ,(j86 
>494,069 
-322,767 
*66,467 
>43.070 
213.221 
313,071 
309.334 
633,364 
723.267 
30». 6*6 
634,161 
700.206 
949,995 
71.073 
92.209 
646.037 
666,144 
336,699 
133.907 
297,336 
770,466 
411,997 
303,974 
353,104 
•120.116 
32.016 
122.213 
326.260 

431078 

-20**4S6 

•260*333 

>386*333 

197.061 

•193.219 

>432*373 

*703*993 

•66.771 

-49,996 

•170.JT1 

•103*363 

139.913 

44*041 

696.234 

136.266 


63MB 

-4.031 
-71,619 
•4.822 
•23i;4l3 
•383,344 
•311,136 
-501, 5$i 
-7*,o7l 
-274,3*9 
••64  7,1»« 
-366,391 
•194,372 
•206,819 
•204.265 
64,45$ 
373,974 
749.472 
614.3i¥ 
776,48* 
464,68* 
161,936 
493  >43 
693,913 
674.72* 
424,87? 
311.321 
310*30? 
184,333 
133,10? 
253.624 
•669,161 
-64*911 
•60,343 
•1199, 34* 
•990,40? 
•1639,423 
—244,063 
127*66? 
139, 2l? 
429,917 
730.46* 
•119.309 
•2.10' 
37 *37* 
673,646 


homo 


327.543 
*11 .693 
379, 8»9 
243,617 
•232,442 
•*31,673 
.266,053 
-*74.612 
-691,763 
•339J766 
•133,749 

•132,767 
•16,049 
.302,244 
-164,709 
110,267 
12.522 
202,799 
401, *»6 
680,049 
614,040 
213,788 
S9*4j9 
409,630 
779,634 
629,370 
419,062 
*08, 0?9 
•o9,*16 
263.097 
llS.Toi 
760,694 
*99,6*7 
511.109 
316.214 
*37,012 
101.3*0 
36.116 


•421,066 

•219,133 

324.119 

-1033,940 

203,173 

•230.392 

350,070 

?26.6*7 


190MB 

1117.362 

702.039 

838,737 

658.946 

274.266 

334.365 

.143,699 

-62.131 

.220.715 

.480.911 

-261,330 

-383.004 

-150.720 

•200,472 

8.319 

130.060 

•136.367 

250.736 

442,094 

300.663 

177.633 

152.409 

-94.864 

•149.976 

40.992 

•206.111 

•190.843 

.ha  all 


67*604 

-48*934 

201*416 

43o*246 

1069*367 

124o*OOl 

1223*007 

1642*629 

1360*172 

1131*915 

839.275 

631.750 

1133*464 

460*057 

93.677 

•369*307 

139*704 

-112*373 
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30M0 

215,378 

*•1,0)0 

196.005 

649,337 

485,883 

«*,3TJ 

227,858 

1172,141 
1181.034 
1*35.023 
1**0,263 
1524,480 
1260,248 
1123,290 
•76,000 
•63.001 
704,061 
*37.020 
104,070 
•150.000 
•108.491 
-22* .016 
-227.347 
64.448 
-26.337 
•0,180 
.174.003 
.150,081 
-♦36,296 
•302.004 
-353.031 
•335,670 
-611.001 
-700.134 
.064,986 
-.867  .940 

-SI?*9*1 

*••1.070 

•1227.743 

•1041,020 

*1270.310 

-083.609 

•1001.347 

.080,778 

-•50.829 

•**2.580 


66MB 

116c  386 
26*:i03 
TSe,o3* 
927.805 
1111.726 
99j,B22 
1130.634 
1104.002 
1256.426 
138i. 06J 

1236.011 

421.100 

573.711 
274.040 
801,134 
73  i. 271 
227.902 
187.736 
239.421 
210.080 
7.643 
-365.936 
-751,32* 
.850.162 
.944,121 

-981.490 

.731,602 
.286,871 
37,561 
.631,606 
-271,752 
-576, 41$ 
-770,482 
.1123.935 
-750.432 
.1420,421 
.621,070 
-754,824 
-416.80$ 
*386,840 
-424.70/ 
-450,830 
-517,50$ 
-500,164 
-655,292 
-456,626 
-227,354 


110MB 
402.863 
438 ,463 
1628,848 
1836,062 
874. 5T1 
15*2.330 
1®13,171 
1653.610 
1$I5,9s0 
1*03,004 
lf09,2io 
1*21,606 
2202,77s 
»««,!« 
1*90, $30 
1J87.843 
$25,361 
*13,869 
*32,560 
220.7/2 
-|19.2o7 
-774,380 

-1*36. 0J3 
-1763,7-1 
.1933,961 
•2108,462 
-1670,1*9 
-1700,3/3 
-2il7  418 
-2207,468 
-2070,262 
•1724.062 
•?lt,10T 
•1005.731 

•902,108 
-706,108 
— $68 ,3o4 
.522  024 
•573,867 
•464,469 
•|74,*O0 
•275,424 
•155,435 
270,407 
Ijl.lJl 
133,441 
♦00,7/3 


190MB 


•174.195 

280.006 

742,403 

1312,361 

769.146 

553,552 

2005,783 

1647,380 

2081,016 

1343.161 

1592.604 

1939.033 

1643,730 

1170,612 

710.467 

632.80* 

517.101 

159.510 

•328.162 

-♦02,301 

•♦73,678 

-350.476 

-531.730 

-1490.726 

-2270.510 

-2727.109 

-3061,310 

.3303,873 

-3572,242 

-2871.008 


-2424.086 
-188) .806 
•125o*Ol3 


-990,654 
-1071,767 
2*1.812 
-375,492 
-969, 454 
-850.055 
-660,366 
•1099,260 


-706.898 

-160.159 

707,680 

574,092 

560.389 

543.517 
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38*8 


•530.398 
•597 .411 
•632* «90 
-1534.196 
•1390.650 
•1075*595 
•663.851 
477.195 
190.664 
47.797 


42.057 

242*369 

52.161 

171.767 

150.934 

207.562 

156.367 

301*587 

162.410 


-38.3' 
7  V  *9  i 
480*61 
250*01 
24l*5< 
•74*35 
•**7.3] 
•86.9« 
<24. 6* 
•252.81 

36*67 


•336.440 

•21**475 

•481.980 

•363.938 

•639.664 

•460.471 

•589.635 

•235.636 

•304.638 

•114.657 

•317.789 

•353.338 

•311.684 

•624.444 


65MB 
•375.301 
•421.611 
345,978 
•115.677 
4.83* 
4«5.484 
9/. 893 
823.02’ 
•636.0*9 
lo.l’i 
-166.794 
10.92? 
*16.074 
395,494 
470.617 
32/. 13* 
379.182 
•165.914 
•95.6*0 
•145.652 
•8J.864 
•455.957 
•585.137 
•47*. o7* 
•251.03* 
521. 16§ 
27^.592 
•734.50* 

•346. li/ 
•311.007 
•3U.02? 
•20o*95u 
“240.18* 
•66.110 
•62.399 
•82*.4flf 
•6,83? 
271,8*7 
-175.981 
8*2.752 
295.018 
779.274 
549.0*7 
265.78 It 
7*1.486 


110Mb 

315.558 

51.329 

•97.250 

♦87.9s6 

225.059 

*43.7/9 

?42.3o5 

268,393 

*38.433 

•77.794 

5.791 
•594.967 
•2*6, 5o3 
-250.469 
•63.528 
249.109 
884.843 
839.413 
61,550 
165.829 
298.5*4 
108.5*4 
•24.636 
•637.2143 
•570.513 
•1065. 9t8 
-629.15* 
•5o3,*u4 
-666,748 
-346,396 
•440.7*6 
-173.510 
-680.541 
•♦29.453 
-4*1,214 
•164,6/1 
-178.133 
♦53.896 
5*1.119 
568.749 
♦73.490 
362.662 

Ml. ♦*! 

761,003 

569.019 


190H8 

613.338 

987.468 

•198.605 

1183.986 

•663.930 

*56.644 

1899.840 

•539.11* 

139.880 

•735.853 

•♦15.112 

•855*665 

-lOSl.417 

-935.446 

•699.494 

144.312 

378.886 

574.958 

*89.166 

744.842 

-17g.862 


141.888 
•109.868 
•748.586 
-582.459 
-687.687 
-815.709 
•348.915 
•603*168 
-382.2*1 
•4*0.995 
•589.746 
•470.834 
-866.553 
-63,427 
819.676 
1*8.817 
380.083 
552.487 
608.997 
79 0.089 
538.1*5 
1063*007 


TABLE  A-3.  U-EASTWARD  WIND,  LATITUDE  24  DEGREES  NORTH 


PRESSURE  LEVEL 


30*0 

1004.992 

«03*813 

1135*405 

ia«f*S33 

1169.189 

13oT.868 
1370.967 
1911.149 
1890*673 
1896.164 
779.073 
613.037 
884.390 
638.046 
386.610 
883.076 
191.606 
-480.619 
•399.884 
-989.896 
•799.843 
•1700.440 
•1880.619 
•1886.974 
•1714.341 
•1980.347 
•1699. $97 
•1691.186 
•986.880 
•1144.847 
•1086*466 
•968.989 


•308*014 

•191.091 

•695.149 


•166.979 

•184.033 

•947.390 

•899.430 

140*808 

-418*814 

•678.440 

•781.030 

•1090*094 

•1077.348 


6988 

1 1 0M7* 

19OM0 

1834.420 

8417,166 

3152.145 

1030.8?! 

8*27.966 

333?. 9i0 

1104.906 

2653. 7«Z 

384?. 958 

1436.066 

2*00.003 

3616*870 

1309.434 

2?40.6&S 

3558*456 

1336.040 

2*47,091 

3955*007 

1174.967 

3039,349 

3092.499 

939. 5?i 

3*30. ?00 

4305.30? 

1480.97? 

2046,295 

4807.051 

1009.87*! 

2060,363 

4646.441 

1006.20“ 

3013,4«9 

4109.010 

1066.740 

3340.938 

470Q.359 

1890.996 

3529.9*2 

4437.730 

1994.019 

3450.3*9 

4640*00* 

1736*99? 

3420.0*0 

4102.043 

796.980 

3003.502 

4837.703 

064.108 

3169,364 

4399.390 

1006.618 

8091.460 

369].066 

006*407 

3140.760 

3866.260 

911.38? 

3133  .?« 

3471.404 

991.996 

3112.298 

3735.254 

616.307 

3403*340 

4019.003 

607.100 

2013.065 

396?. 071 

400.607 

2404,244 

3994.120 

086.804 

2057,464 

4133.294 

603*973 

3107.003 

4047.441 

134.004 

3017.686 

4310.312 

104.8?* 

3109.036 

4207.965 

209*514 

3414,677 

4495.344 

-.158 

3313,030 

4431.121 

40J.411 

3020.069 

4633.375 

414.301 

3*73,5 J4 

4514.324 

784.166 

3001,931 

4009.144 

676*406 

3644.761 

518?. 020 

081.676 

2727.4*4 

4345.355 

1064.101 

3904,0*7 

5011,031 

1766. 3?6 

3278.461 

4531.199 

1014*706 

3010.308 

4498.010 

710.00* 

2070.739 

4410*420 

936.194 

2?30.094 

4908.941 

434.800 

1049,369 

4152.477 

271.05* 

2131.050 

4504.240 

036.467 

1940.664 

4509.133 

493.03* 

8098.319 

4218.465 

747.30; 

8479.741 

3004*131 

634.024 

8479.260 

4277.480 
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TABLE  A-3 .  CONTINUED 


3»l»8 


-10B6.ee? 
-1311.040 
-1345,907 
-1*37.267 
-1381.714 
•1820 .328 
-1*1*. 626 
•19*9.300 
•1043.627 
-2223.054 
•1827.631 
•1086.903 
•1080-971 
-2075.076 
•1081.019 
•2997.920 
•2792.719 
•3037.968 
-3120.718 
•332* -9 14 
•332*. 970 

•3036.291 

•2988-394 

*2*23.987 

•2*39.07* 

•2*10.353 

•2890.299 

•2997*384 

-2824.131 

•2692.098 

•2727.931 

•2697.407 

•3210.673 

•2089.991 

•3193,098 

•3292.199 

•2786.690 

'3026.691 

•2886.697 


•2729.189 

•2647.101 

•2*20.390 


-2*66.39* 

-2286.02* 

•2001.679 

•1886.096 

-1989.071 


6988 
60*. 201 
26*1*90 
274,6*3 
90.28* 

*o*;oi* 

0,827 

187.893 

•*6*,326 

•555*116 

•670.505 

•516.320 

•624,11 5 
•824,281 
-106*4789 
-13*64*92 
•1931.280 
•1513.951 
•1537.50V 
*1707,32® 
-18624334 
-233*4*04 

•2330. *2® 
-26**. 37? 
•2*7*. 63« 
•2402.504 
•1007.631 
•U02480* 
-711.55* 
•1017.556 
•1*06.017 

•13*«.Ol0 

•1*53.16$ 
•03*»55« 
•iii*»«oJ 
•103**162 
•876.103 
•087,785 
•03«*9l2 
•715. *64 
•*04,80* 
-10*,352 
•583.566 
•204.503 
•1 0. .52* 
3*. 827 
477.384 
*36.12® 


1108* 

2373.051 

24*6 ,499 
2*35,112 
2076,744 
1*0*. 049 
1150,000 
1314,617 
1302.211 
1592.3*8 
1230.547 
S00,B*8 
•16.70T 
276.118 
35.567 
•428.8*4 
•9,5*0 
•*6.7*7 
-287,4*6 
-588.0*8 
-4*3.170 
•418.279 
•7o6,563 
•619.667 
-1123.209 
-1319,181 
-076,55? 
-§35.095 
•?62,0l8 
•137,004 
107,715 
132.010 
75.0*6 
307.726 
•00. 0*2 
♦O7.0*o 
1504.002 
0*0,15* 
•01.826 
1 150.2*8 
011.0*7 
1105,221 
1028.03* 
1471,232 
1022,4*1 
1964.766 
l25l,0»7 
1408.243 


IOOHO 
30*6.201 
3376.450 
3065.090 
2625.135 
2051,600 
1279.575 
1296,8*7 
1581.650 
1511. TOO 
1410.325 
1335.5*0 
078.085 
531.255 
0*0.453 
027.761 
736.347 
082.937 
9*2.434 
•  u.105 
488.614 
869.630 
667.006 
268.302 
•101.006 
53.442 
252.718 
530.5*1 
578. **• 
217.871 
579*036 
65}. *57 
1222*867 
1831.7*7 
1734.585 
1053*473 
1716*2*5 
1361.9*8 
1820.810 
18*8.618 
1800.603 
10*9,068 
1718.785 
1050,838 
1101.738 
1680.800 
1767.1*8 
10*0.671 
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TABLE  A- 3 .  CONTINUED 


38*10 
•14*8.401 
•1 043.9*1 
•▼*4.5ji 

•54«.o*' 
•471  .*4] 
•2l9.*0] 
*4**331 
•333.843 
1»1«4!< 
145.041 
44*.  1*1 
4*2.031 
•97.791 
4?4.404 
«8«.iio 
510*800 
*45.298 
768.040 
**7.372 
•14.930 
534.639 
250.050 
•54.801 
1091.715 
120*.950 
1147.067 
•08.790 
493.636 
1040.208 
•*0.o*7 
•81.614 
•*5.760 
▼•*.7*9 
1294.00] 
775.794 
1050.087 
•*7.743 

•o«*6is 

323.1*9 

870*803 
•1.9*0 
215.974 
2*7.752 
112.812 
4*4. m 


•HfcSSuUfc  LEVEL 


**M8 
227.44V 
524.230 
807.479 
77*.43» 
1210.763 
96«, 0*4 
1706.61* 
•16.614 
868, 8lS 
344.39* 

7*5.70* 

711.941 

1126,41? 
11*0*00* 
777. *97 
•74*008 
10*2 .*•! 
•35.501 
*40.34* 
77*.84t 
910.730 

147*. **T 

*93.074 

450,94* 

2M.*97 

510.984 
1*71.354 
l2«kf 912 

78* .458 
1031,9*7 

880.984 
•**.23| 
*88.578 
1*7,774 
4*0.14* 
434,6*9 
171.77* 
8*4.317 
45o»2*& 
09y.*lt 
*43.87* 
934,627 

1212.5*4 

1134.411 

157o.22i 


11  OMb 
1268,442 
1395.087 
1ST*. 7«i 
1277,23* 
1470.130 
2653.31* 
2106.4*4 
1*87. lil 
0140.269 

2190.000 

2170.538 

2?]0.0V* 
2035, *33 
2215,58* 
2147,0*6 
2332.413 
lf73,9»4 
2396,469 
2382.8k* 

l8l0.«3* 

1?57,4*9 

1*77.5*1 

2118.921 

2245,6*7 

0546,499 

0273.055 

1086.5*0 

2305,000 

2220,197 

2^13.3*7 

1740, *05 

1626,969 

1645,233 

1*72,8*5 

1 396, 5/8 

1*88.7*6 

1321.837 

1832.0/6 

1*24.2*5 

2011.4/0 

2318,023 

2074.418 

2*35.473 

2*33.474 

3173.528 


190MB 

684.915 

1503*412 

1256.561 

2252.164 

1803.1*3 

2309.148 

0350.089 
3057.7*3 
2320.4*1 
2685.757 
2546,930 
3115.416 
3330*645 
3579.07* 
3365*855 
3308*01* 
3354. 0*8 
2676.079 
2482,479 
0945*897 
2731.574 
2221.1*7 
1*96.0*9 
2170.740 
0350.817 
2599.513 
1948.804 
2070.05* 
1741.511 
1*39.479 
1706.1*8 
1657.039 
1*05.022 
1720.999 
1**5.262 
1746.04* 
1817.08* 
19*4; 72i 
18*9.14* 
205*. 712 
22 *5.712 
2449.041 
2929.778 
2849.670 
278].898 
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TABLE  A-4,  w-V£Rt;  • 


PE  48  DEGREES  NORTH 


39*6 
Mil 
»««{ 
1.699 
10)1 
1.264 
MM 
I  *M9 
MIT 
3*179 
♦M 
•  e?? 

.$79 

•Uo« 

•JT| 

2*763 

-.2*3 

•1.916 

•973 

.*94 

.995 

1. *17 
*210 

-.999 
1.9W 
1*699 
1.949 
♦  •416 
3.9*1 
2.794 
♦•119 
9.960 
9*001 

2. (|7o 
1  .§96 
1.129 
1.369 
2.390 
-.393 
-.369 
2.719 
1.972 
-.261 
•*292 

.960 

1.437 

-•0*9 

•979 

*619 

1»€»7 


4**N» 

.92* 

•  94» 
4.39* 

•1.99* 

«91» 

«.*♦* 

2*2*4 

1.77* 

>1.4*W 

•  4ll 
•  •3*9 
2.8*3 
2.621 
2.319 
— 8*« 
i.07l 

.132 

1.223 

1*29* 

1.40$ 
1.200 
2  .4*S 
.97* 
1.47$ 
*.123 

k. 21? 
».90« 

l. 129 


.3*1 
•  1*0 
.♦*# 

•  oio 
-.0*1 
'it  i* 
1.9-0 
-.1*0 

-1»*44 

•  7  *9 

•  4«* 

.♦•3 

.6*9 

-.020 

•2.2)) 

-.•ll 

.*22 

-1.092 

•1.3|0 

Ji2 

-.Oto 

•  1)0 
1.679 

.696 
•  3*2 
.0/9 
2.113 
1.0*6 
3.6)1 
♦»2o4 
1.749 
•»6<0 
1.013 
1.6*1 
1.799 
1.1*7 
.270 
1.632 
.6)0 
-.022 
.0*3 
.0*7 
1.3*3 
.*?B 
.*•2 
l.$J* 
.0*4 
•1.613 
— 1 .3*1 


190MB 

.730 

••300 

-.139 

-.411 

•1.0*7 

>.470 

•  5*0 
1*649 
1*200 

.349 

.372 

•  017 

-1.3*8 

-1.472 

•  12» 
-.107 

-1.932 

.128 

1.1*2 

-1.463 

-1.078 

-1.055 

-.737 

.541 

*006 

.040 

—on 

•  70S 

•  756 
2.61* 
3.186 
2*320 

.7*7 

1.353 

1*081 

1.466 

1.731 

.762 

-.246 

.446 

1.460 

1.880 

•  673 
.321 

•  016 
.912 

1*703 

1*986 

2*10* 
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TABLE  A-4.  CONTINUED 


38VR 

.719 

•2*57? 

-2.55J 

1*055 

»4?3 

•2.125 

-l.BlT 

*•830 

-.031 

-.*91 

•1«1 

•1.038 

•1*310 

*•435 

-3*531 

•3.7*0 

•2.74* 
-1.89? 
•2*3*9 
•4.774 
•3.030 
•1*9*0 
•4.474 
-3.0*4 
-1*833 
••422 
•1.247 
•2.01* 
•2.447 
•2.410 
•1*930 
•3.19* 
•3*342 
•9*050 
-3. f*o 
— «do 

••114 

•1*295 

•1.4*7 

•1*311 

1.0*9 

1.082 

2*233 

2.100 

1*447 

1*381 

2.040 

2*000 

1.737 

1.237 


OS  89 
.«H 
2.7  91 
•••♦*$ 
-1.50* 
2.542 
*02w 
'••485 
■••99ii 
«.9B* 
-l.Olt 
.52* 
•2.519 
*3.224 
•  24§ 
•3*571 
*9.560 
-.0*7 
•1.680 
•1.951 
*3.409 
•4*259 
•3.51* 
•4.41* 

-.214 

.Oli 

•1.01* 

•2.905 

•2.664 

■*022 

.505 

•1.702 

•3.372 

•.i7i 

•3*100 

'*•070 

*3*049 
•  009 
-.740 
-.94* 
.450 
2.024 
2.902 
2.«1* 
3.35* 
2.004 
3.5*8 
4.009 
3.135 
.721 
2*090 


110*8 

1.019 

5.390 

3.035 

•1.581 

-.242 

•  8v? 
.724 

-.593 

•1.110 

•1.197 

-.679 

•2. 3vo 

•4.152 

-2.1*1 

•.790 

•l.*97 

•2.598 

•3.113 

•1.485 

-.342 

•2.734 

•2.032 

•.9»5 

•  5oo 

•1.040 

•2.4*1 

-.1*6 

1.5>o 

1.4l* 

-.4*3 

•I.«*l 

*3.0*0 

•2.6*0 

•1.195 

•2.203 

•3.0o* 

•2.404 

•*?91 

«5iO 

’I'3 

.<47 

1.582 

2.2*0 

2.8*3 

2.4*7 

2.8*5 

1.502 

1.400 

.705 

l.?oa 


190M8 

-.039 

•j.899 

•1.727 

•  .720 
*049 

-.749 
•1 .067 
-.905 

•  •374 

•  »089 

•1.617 

-.910 

•  224 
*j«444 
*2*48q 
-1*224 
•3*245 
•3*840 
•2.548 

1*002 

.503 

•  342 
1.774 
-.506 

•1.539 
-.930 
-1.008 
•1.315 
1.155 
1.104 
-.989 
-2.259 
•4*050 
•1 .986 
•»305 
-1.218 
•1*027 
•1*791 
-.547 
-.105 
—  .564 
-.137 
.423 
1.125 
1.374 
1.770 
.027 
.782 
.200 
.059 
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TABLE 


A- 5 .  V-NORTHWARD  WIND,  LATITUDE  48  DEGREES  NORTH 
PRESSURE  LEVEL 


38f*B 

-486.752 
-613.(83 
-343.640 
-704.730 
-491.703 
•141.254 
-557. Jo7 
•36»o63 
301.639 
209*690 
867.365 

724.355 
501*260 
449.073 

1379.262 

14*8.560 

1064.1*6 

961.462 
478.2,2 

451.066 

8o.l?6 

-269.4*5 

-366.325 

38-363 

-102*392 

•104.136 

-125.860 

-59.508 

-111.079 

-255.376 

353*263 

67l.l*9 

850.463 
1053*069 
1391.323 
1179,369 
1673,922 
1887**63 
1860.992 
2089.67] 
2197.867 
2383.460 
2093.0 65 
2229. *44 

1716.355 
1434.292 

952.145 

649.276 


69HB 

110*6 

190M8 

-83 / .269 

-712,526 

-365.742 

•874*426 

-193.707 

•61.441 

•544.487 

-403.736 

•261.201 

-413*764 

-145.310 

321*135 

-476.844 

-230. 606 

2*5*078 

-61V. 666 

•523.562 

-416*6*7 

-640.543 

-*04.210 

-148]. 890 

-6.477 

-771,468 

-1620*387 

308.057 

•?o 1.809 

-844.254 

—  1  Jo .513 

-4)6.404 

•361.969 

*183.427 

143.469 

•45*113 

-178.153 

394.569 

213*3*7 

-294.794 

772,404 

469.3*6 

748.254 

677.145 

513.153 

1*64.761 

397.1*6 

1176.887 

596.829 

687.935 

848*087 

346.44* 

1161.358 

444*848 

518.144 

774.142 

923.960 

1031.48* 

638. ?04 

821.307 

593.604 

519.524 

1113*553 

414.513 

460.374 

1105*21* 

-133.149 

652.3*6 

696*08] 

-17*. 80* 

357.273 

69] .90 7 

-6/. 598 

315.173 

753.271 

-26.723 

‘59.921 

179*460 

•486.004 

-635,*86 

•557.694 

-1054.956 

-ll»7,90S 

•1516*417 

-1033.493 

•  IU6.48J 

•1785.645 

•1328.565 

■-1**3.264 

•2660*673 

-1176.734 

-1*42.740 

-3078*460 

-476.387 

-1*67,0®* 

-2876*8*6 

-35**72? 

-1571. 5o2 

-1719.837 

-530*78] 

-1064.823 

•984*756 

37.956 

-699,806 

-979*056 

-124.76' 

•305.487 

-611*890 

680.576 

130.373 

-424*251 

1280.971 

249,970 

•94*076 

1353.37* 

7*6,403 

208*912 

1677.170 

1205,149 

598.889 

2288.720 

l702.496 

619*348 

2594,13] 

l6o8,974 

704*371 

1910*204 

1248.962 

1220*798 

1836*873 

1656, *66 

1697*694 

1340,053 

1**2.023 

1500*282 

1023*356 

l650.1*4 

1500*337 

552.31I 

1342.116 

1101*53* 

69*|  ,40* 

8*5,474 

6*1*105 

1143*328 

582.170 

-88*541 

1348*460 

390*560 

-704**50 
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TABLE  A- 5 .  CONTINUED 


36MB 


— v— 

697. 756 
670.744 
♦52.074 
3o«*294 
112*644 
•137.li* 
•*46.277 
•7*3.742 

•1344.035 

•1572.843 

•1518.249 

•1*47.635 

•2150.495 

•25l7.4?9 

•2590.061 

-2791.149 

•2*04.138 

•2**2.934 

•2772.2*1 

•25*6.*92 

•2349.150 

-1987.846 


-1*28.461 
•1*21.21* 
•130«* 946 
•729.3*1 
-315.60* 
-130*678 
571.(42 
775.1*1 
1*30.146 
15*3.210 


1551.650 

1*78.863 

13*8.911 

1201*544 

1296.712 

1027.950 

496*394 


•137.563 

-♦69.9*1 

•669*639 

•685.788 

•718,130 

•884.595 

-876.221 

•1081.757 

-6ll.o*9 

•502*657 


65M0 

UlO.68* 
6*0  *0** 
534.921 
*7*. 7*2 
274.70* 
•40.462 
3*6.731 
•*0,7*3 
•261.0*2 
•57/. 196 
•956. 9Bw 
•1206.133 
•lias. 57? 
•1730.14* 
•1693.20* 
•1664.152 
•20»3«*7l 
•1670.3*3 
•1690. 555 
•1546.13* 
•1176. *0? 
•lo9*«2*7 
•1276.550 
•1336*335 
•526. *71 
•19.663 
•106.22? 
26/.29| 
229,6*5 
6*1. 6*7 
1024.512 
1156.30* 
1240,6*3 
1595*029 
15*7.385 
15*3.690 
I2*6.*l7 
*17.1*0 
213.998 
-154.527 
•♦34.276 
•641.67* 
•1020.22? 
•1023*13* 
•953.935 
-1039*166 
•733.623 
-574.593 
•076.5*5 
•79*,253 


110*8 
•107.197 
•75.691 
•216.799 
•22*. •*? 
160.278 
160.632 
32*?24 
“131. lo° 
-261,676 
•356.218 
•*27,145 
•327,2/0 

•21*. 536 
•270. *04 
•*11.2*9 
•*94.6i7 
•552. *49 
-*73,942 
•525.484 
-*55.73! 
**21 . 0«6 
-«75.8i9 
•*07»*wi 
•55.0«7 
-311.597 
•♦10.145 
-*19.19* 
•5o3.*86 
•i*7  «5|9 
*63. *92 
12*3,700 
2**7.263 
2328.672 
2103.535 
1682.959 
1*21.011 
1**6,115 
1137,572 

•*n.*oi 

•1**1.028 

•1*01,546 

•1*09,060 

•1150,064 

•1764,206 

•1358,9«3 
•1*67,429 
•10*0,127 
•1*31. *12 


190MB 


-706.750 
•1*0*9*? 
386.71* 
399.611 
*59.5*9 
-36.467 
•274.99* 
•71.712 
•2*4.13* 
•13o.052 
-1*3.271 
219.192 
201*106 
365**73 
517.909 
887.009 
841,069 
374.01? 
**4*552 
-55,,. 772 
•755.0*6 

•*39.697 

•*75.2*5 


-346.730 
235.374 
443*49* 
327.0*1 
•292*209 
•*05»l9* 
08.313 
1235.006 
2*81.412 
2994i 14* 


2*42*665 
2017.066 
15U.583 
1634.402 
1*03 •*!* 
*29.177 
131.206 
•844.442 
•1701,84* 
-2*93. *25 
-2668.814 
-2347.2*? 
-2232*905 
•1800.050 
-1283.947 
•621. 012 
•**6»*12 
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TABLE  A-6 .  U-EASTWARD  WIND,  LATITUDE  48  DEGREES  NORTH 

PRESSURE  LEVEL 


38KB 

65MB 

110H8 

190MB 

4191.091 

3698.080 

2202.276 

1103.367 

4390.5TB 

3613.0*1 

2*69,048 

1203*16* 

4610*937 

3497 • 0** 

2*07.782 

1*44.852 

4676.773 

3624.03? 

2*90.0/7 

1703*1** 

4567.973 

373/. 02* 

2*40.2/3 

2138.*** 

4102*900 

372(j.o3? 

3046.7/9 

2808*137 

4464*246 

350/. 01’ 

3177.4/9 

30*4*6*0 

4539.393 

336o.32& 

3057.720 

2*11*1*7 

3829.2,2 

3960.15? 

2*10.028 

20*2»**0 

4444.664 

4162.36/ 

2506. 872 

22*0*9*1 

4318.840 

3904.288 

2716,216 

22*7.*** 

4479.746 

3969.674 

2682.899 

25*2*41* 

9343.664 

3592. 0*5 

3065.1*6 

2819.876 

3192. «4T 

3l44.Q6t 

3*98, Buq 

30*3.957 

4666.723 

4063.04* 

3563.908 

31*5.431 

4541.629 

4796.605 

3342.363 

3883,189 

3328.890 

3937.093 

3647,417 

35U.240 

5664.203 

4038.043 

4127.496 

3817.920 

3796.195 

3970.90* 

3*2*. *11 

3214.307 

3792.359 

434.. 0S» 

3?8 1.646 

3273,964 

3438.625 

4244.07* 

3371.201 

3210*469 

9389.297 

4306.410 

3373.610 

2942.233 

3300.789 

3817.772 

3188. 09B 

2336.04* 

9210.258 

3968. 88? 

3166.6*6 

2164.075 

3066.730 

3998.135 

2837.263 

1784.018 

5430.949 

439o.937 

2j|3.7uo 

1673*802 

5178.340 

4230.140 

27  02.632 

1883.186 

9194.062 

403/. 356 

2*33.9*6 

1807.863 

4970.293 

4004.9U 

29o9,8jB 

20*0*30* 

S2o6.766 

3399.062 

3037,992 

2224*030 

4597.023 

3337.82? 

2*92.1*0 

1847.312 

3969.970 

3786. 32U 

2*25.573 

1379**70 

4489.582 

3683.422 

2*59.649 

1465.031 

4370*273 

3609.38* 

2328,032 

1647*684 

4743.399 

3841.53 i 

2*1*. **l 

1873*793 

4862.469 

3608.131 

2?92.738 

1839.949 

4762.070 

*132.250 

3078,266 

2287*356 

4648.430 

*614.043 

3532.143 

2780*634 

3017.059 

460/. 62? 

3668,136 

336J  ',333 

3084.031 

4949.30? 

4523,981 

3877. *23 

4453.380 

538,1.973 

5162.3*3 

394],760 

4487.566 

3801.105 

5339.7*2 

4514.333 

4423.422 

5023.340 

5369,622 

4901,078 

4221.631 

5996.590 

3785.9*2 

3800*480 

4089.397 

3706.512 

5*80.539 

6129.71* 

3900.793 

5276.00? 

6431.5*4 

6747*335 

4003*317 

4716*32? 

6379.219 

6874*363 

4094.741 

4831*703 

6*83,973 

70*3*2*7 

4031*332 

4008*21? 

6266,018 

6199,004 

A-15 


TABLE  A- 6 .  CONTINUED 


38mS 

3903,6*7 
3583.167 
*026,595 
*069. *?5 
*0*6,1?* 
3877,863 
**71,516 
*187, 623 
*175,51? 
♦203,836 
39*1,530 
399*, 721 
*059,769 
*035  e«3 
*316.137 

♦602,3*0 
5090,0*1 
5613. ?15 
5674.488 
5901.053 
6138.465 
6269.5*7 
6273.996 
6*93.336 
6201,666 
5593. 1?9 
55?7.*33 
5532.16* 
5*57.203 
5160.367 
5275.236 
*115.332 
5523.785 
5*72.312 
55ol. *22 
5755,8*8 
5269,016 
5148.073 
5283,162 
5519,73* 
5*71.361 
515* ,539 
♦896,16* 

♦884,520 

*630.812 

*729.730 

*3*1,16* 

*172,878 

3928.0*2 

♦*30,875 


6SMB 

5529  00“ 
527«;95] 
*66«;937 
*619.082 
*88«.U? 
**20.736 
*2** .25“ 
*336,91« 

*369,691 
*281.10“ 
*396 ,672 
*174,73* 
*34*. 00* 
*S5j*27i 
*378,55? 
*73*. 34* 
*698,387 
*866,437 
50U.781 
♦966.12“ 
5006.133 
532* , 34* 
528o, *6“ 
*7U,25| 
♦2*6.27? 
♦566. 5t“ 
*679,996 
*521,25* 
*661,590 
*326,35? 
*505,262 
*656,281 
*354.58“ 
*316.  -iO» 
*276,786 
**8i;o3l 
*610.305 
*581,551 
*198*781 
3931  562 
3808.622 
3839,253 
3789,275 
3738,026 

3324.256 

3721.55* 

3289,562 

3851,966 

3730.461 

3638.02J 


110MB 
6lOB  6u2 
5671*148 
5*51,215 
*771.562 
*f3i ,260 
♦*77,089 
*533,690 
*088*080 
*686,26* 
3608,4*6 
3?37 ,2si 
3575.9*8 
3|85,2>7 
3f*0, ill 
♦^97*707 
*076.?oi 
♦437.68? 
*506,066 
*801, lu9 
*«2*,83Z 
*820, 7 ** 
*?39,6*5 
*752,605 

5609, 39l 
56*9,9*) 
♦J01,695 

♦?03  2s* 
3694,58* 
3572,029 
3338.372 
3500,782 
3*03.682 
4137,8** 
*198,309 
3500,0*0 
3292,7j* 
3ll0,6l7 
3233,190 
3239*317 
3370  812 
3*38.712 
3525,799 
3270.1*7 
3151,890 
3027,16* 
3082,3*8 
2*69,730 
3088,4*3 
2475,1*3 
2*63, 7*6 


190MB 
5370.070 
*875*629 
*866,668 
♦7*3,092 
♦58 1 ,887 
*183.168 
3970,462 
3553,830 
3*03,270 
318*. 876 
3067,562 
3005,989 
263«,619 
2*20,193 
3030.052 
30U.633 
3750,538 
*215.195 
50*1.281 
5891.352 
5*9?, 7*4 
5760.031 
5160.766 
♦690.672 
*170,32* 
3570.256 
3587,857 
3661,751 
36*6,537 
2990,965 
3219,4*4 
3*8*. *56 
37**,*03 
3**4, *0? 
2927,968 
2615.622 
2199,271 
2374,030 
2902; 039 
3*08,979 
3580.606 
356*,823 
3283,421 
278*. 330 
2353,686 
2*65. 16* 
2287; 169 

2201.856 

1796;338 

16*0.466 


.  rl’3 
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A-7.  W-VERTICAL  WIND,  LATITUDE  72  DEGREES  NORTH 
PRESSURE  LEVEL 


30*0 

68*8 

-.453 

•  •460 

•360 

1*281 

-.151 

.022 

1*208 

•  443 

•  098 

.77? 

1*031 

•  70* 

•  941 

1.314 

.550 

•  12f 

•510 

*124 

•  099 

l*09u 

1*329 

1.16* 

1.073 

1.07< 

1.150 

1.2?7 

•411 

.981 

•  106 

,27t> 

•  160 

•  182 

.368 

.0** 

•323 

.109 

.603 

.23! 

.777 

•  57» 

1.414 

.6*0 

•39] 

.782 

•179 

••lit 

-•261 

•99* 

••241 

•  16? 

•971 

•*•993 

1.  57* 

•«3i 

••334 

••031 

1*159 

1.202 

•4*0 

1.3*2 

.489 

•*1  ] 

•  0T9 

••01* 

••14? 

».o«5 

*•!#* 

•  IT* 

••*77 

•••sou 

••6*2 

••01; 

••799 

•*22i 

-.739 

*.96f 

•  •579 

••oi3 

•303 

•  •616 

••369 

.072 

•Mil 

•2.399 

•1.434 

•  •027 

•1*939 

•i.Tti 

••949 

•1.889 

11 0*0 

•  1?1 

•  2?l 
.0'6 

•  *!• 

•  2»3 
■  l«3 
*010 

•  0*6 
•M 
.Tas 

•I*2 

•*oo 

.•*o 

,2is 

•  ois 
•Aid 
,2*2 

••♦It 

-.01? 

•h( 

•  9*2 
••108 
*.T»3 

.1*4 

•  DOT 
.9®7 

•Ji* 

•In 

••o»s 

••1*3 

•  MO 
••0*0 

••1*8 

•  •8*9 
••131 

•  ,930 
*1.000 

••♦19 

•  •*30 

•  •Til 

•  i**t 


190*8 
•  •364 
—  6Tl 

•  773 

•  022 

•  438 

•  0*1 
.389 
•  00* 
»14| 
.813 

.140 

.380 

.023 

.062 

.414 

.292 

.107 

•90S 

.327 

.19* 

.013 

-.360 

•  iTO 
•31? 
.001 
.063 

-.271 

-.167 

-.998 

••400 

—302 

-.327 

-.490 

-.202 

.171 

-.293 

•  •043 
..392 

•  TOO 
-.109 
-.912 

.174 

•1.093 

-.944 

•  •237 
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TABLE  A-8 .  V-NORTHWARD  WIND,  LATITUDE  72  DEGREES  NORTH 

PRESSURE  LEVEL 


38*8 


120. 784 
142,2*2 
.206.976 
.592,119 
.434,6*5 
-635.294 
.169.937 
-44,779 
108.012 
48,831 
.104,118 
.43,896 


.33,282 

.130,880 

-204,307 

244,138 

379,420 

373,583 

311,194 

321,879 

.83,599 

.331.752 

-396,090 

.445,598 

.631,854 

-987.776 

.9o»,667 

.1140.563 

.1291,974 


-1027.348 

.968,870 

•602.440 

-280.745 

.20.481 

732,763 

946.677 

1193,257 

1517,542 

1402.709 

979.511 

604,809 

1087,894 

720,860 


602.349 


485,991 


6bMB 

110«8 

190MB 

18<  638 

151,807 

40.631 

-18 

.300.901 

.344.986 

.435,938 

.261.054 

.195,428 

.788,969 

.670,597 

.693,366 

.424,622 

-632,197 

.923,878 

.650,091 

.563,199 

.799.626 

.858,061 

•?25,?  (3 

.578,424 

28,421 

98,104 

.110.185 

269,901 

111.182 

•198,100 

.174,610 

-473,987 

.579.524 

.338,861 

.623,275 

.741.045 

.328,919 

.738,890 

.778,713 

.542,402 

.1004,399 

•1027.062 

.436,984 

.618,601 

-687.767 

-276,193 

•68,647 

.86,589 

,089 

97,225 

275.770 

64.171 

27.535 

244.555 

-If ,676 

•92,600 

13.850 

-226.892 

.139,017 

.143.521 

-142,557 

129,641 

48.652 

-149,409 

60.818 

290.796 

-81.893 

.144,396 

.77.500 

.76,939 

.105,445 

.468.333 

.204.421 

.356,749 

.469.752 

-274,754 

376,732 

225.244 

-422I92S 

-83;49i 

251,922 

.828,261 

.693  *04 

.519.455 

.848 .26 j 

-§50:?09 

.650,605 

.913,352 

.406,265 

.917.322 

.1144,961 

.959,005 

.908.625 

-974.789 

.769,685 

.736.024 

•381,654 

-334,215 

-487.160 

-100.931 

.164,599 

.351.266 

186.466 

291,799 

90,666 

761,399 

1000,897 

1089.244 

1058,593 

«37:5w0 

758.622 

1 16  f  281 

6#9:i4l 

651,674 

1720,265 

1461,143 

1449.159 

1340,70* 

1262:611 

818.953 

576:923 

943,000 

1410.266 

731,342 

1439,992 

1095. 116 

871,072 

460,193 

574.391 

807,823 

600,998 

1314.257 

S2«, 058 

§59,459 

1147,267 

469.269 

619,675 

1107.010 
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TABLE  A- 9 .  U-EASTWARD  WIND,  LATITUDE  72  DEGREES  NORTH 

PRESSURE  LEVEL 


aiffl 

16o9.4?S 

1347.832 

1147*000 

1079.320 

1230*702 

10««*2O3 

1142*079 

1130*09; 

1421.720 

1974*022 

1724*490 

1633.139 

1914.207 

1402*449 

1302.90* 

1390*992 

1909.019 

1990.903 

2140*300 

2770.790 

3097.969 

3040.12* 

3110.072 

3391.002 

3102.013 

3499.077 

3440*339 

3193.091 

2709.401 

2064.1 69 

1919*090 

1137.942 

040.417 

992.234 

730*222 

994.4^6 

1242.230 

1400.276 

1261.499 

1690.964 

1479.062 

19o6*9*1 

1449.292 

1212*061 

1203*672 


09HB 

110*9 

190MB 

701.073 

*58,946 

-152**75 

741.140 

104.106 

131.967 

610.660 

64. Too 

•440*376 

313*077 

-452.921 

•696*422 

324*96 l 

•434.547 

-716.77# 

120.379 

-462.363 

■814*9q6 

193.63$ 

-3*3.705 

•672*742 

3211*610 

-122.602 

•416*574 

761.624 

388,278 

100*7*6 

978.450 

7 13. 1*6 

414*205 

1179,930 

1030*2*5 

814*607 

1304*969 

1137, ?22 

#21*439 

1330.609 

*11.099 

39). 758 

1390.016 

775.201 

181*684 

1179.555 

429,964 

->84*752 

720.031 

"66, 6l* 

•421*633 

674.019 

350*398 

2*256 

690.160 

£46,606 

51*702 

1043.232 

226.199 

•34*825 

lQ4o«0>* 

225.612 

•56*270 

994 .66$ 

624.066 

204*684 

1220*931 

3l5.?44 

-42*711 

1327*332 

663,302 

8U.T39 

1699.296 

697,677 

210*350 

1497.204 

593,668 

70*692 

1599.055 

630.381 

-316*03* 

1617.970 

539,301 

5*633 

1967.011 

561,000 

-42*511 

1727.510 

*09,234 

360*391 

1460.9*6 

*56,719 

461.486 

1240.94* 

*53,345 

635*3*6 

667.251 

455.241 

374*0*6 

610.077 

676,244 

565*017 

410.167 

322. *«0 

516*376 

419.653 

259.747 

272*01* 

609.556 

672,028 

479.163 

694.500 

2si.6oo 

48.526 

1154.916 

•71,673 

838*561 

919.259 

623,938 

469*355 

1040.021 

644.022 

623*612 

909.6*3 

383.200 

•112*071 

6li.70« 

460.677 

•46*320 

663*054 

-S6,ioi 

•757*967 

4*9.274 

48.304 

-227*642 

6S|.9i3 

116.791 

•158*751 
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TABLE  A- 1 0 .  COMPUTER  PROGRAM  FOR  CALCULATION  OF  a-STATISTICS 
FROM  3-0  GCM  WIND  DATA. 


PRIM  2*»00 

print  6oo 

PBTNT  200,  («UM*4e<J)  ,J*1,4)  »<AV/GW4tt(K)  ,K*1,4) 

PRINT 
PRTNT  600 

PRINT  200,  («UMWP4e<J)  ,J*1,4) «  <AVG*P4b<J>  ,J*l»<n 
PRTNT  77OO 
PRTrT  6P0 

PRTNT  200*  ( (ALPl4e  *J*1  *4)  ,  (ALP248  <  1  ,K )  ,K*  1  ,  4  J  ,  I  *  1 , 99 ) 

PRINT  3*00 
PRImT  600 

PRINT  ?00,  (CNTA148  (J)  *  J*lt4)  ♦  (CM  A<f48  (R ) ,Kx J ,4 ) 

PRINT  3^00 
PRIM  600 

PRINT  20  0, («UMAl4g  < J) ,J»1»4)  * (SUMA24B ( K ) »K*1 ,4) 

PRINT  y?00 
PRINT  600 

PRINT  200  «  (  AVGA146  (J )  »<J*  1  »4)  * (AVGA24H  «R)  »K*I  *4  ) 

PRTNT  4100 
PRINT  600 

PRINT  200, <«CEVAl2 ( J) ,  J*1 ,4) * (S0EVA22 (KT »K*1,4» 

PRTNT  4”*00 
PRINT  600 

PRINT  200, ( (CEVA148 ( I , J) , J«l»4) » (0fcvA24O ( I ,M »** 1 *4 >  » 1*1 *9* > 
PRINT  a«00 
PPTnT  600 

PRINT  200, (A VGA  12 (J) ^v*l,4) , < AVGA22 (R ) ,R*1 ,4 ) 

PRINT  47oo 
PRINT  600 

PRINT  20  0,  (BGl4bC„>  »J*1*4) ,  (KG248  (R  >  »K*l»4) 

PRTNT  7^00 
PRIrT  600 

PRTNT  200,  <  <U72  (I  ,«J)  »J*1  ,4)  »  (DEVU72  U  ,R)  *K*1 ,4)  *1*1  ,45) 
PRINT  24Q0 
PRTNT  600 

PRTNT  200,  (  (V72(I,J)  *J*1*4>  »  (DE  VV7*  1 1  ,l\)  ,K*1,4)  ,1*1,45) 
PRINT  2*00 
PRINT  600 

PR  t  Nl  200,  f(N72(I,J)  ,J«1,4)  »<0EVW74U,n>  ,K*1,4)  ,1*1,45) 
PRTNT  ?«oo 
PRINT  600 

PRTNT  200,  (SU*u72 < J) , J«l»4) » (AVGU72TK) ,K*1 »4) 

CONTINUED 
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TABLE  10 


CONTINUED 


PR T K  T  1OO0 
PRIM  «C(i 

PRINT  ?OCf  fSUMV72 ( J) * J*1 *4) * (AVGV U  ( K » *K  =  1 »4) 

PRTM  ->100 
PRTtJ  A00 

PRIM  200»  f  ?UM**72  ( J )  * J* 1 *4)  *  (AVGW/^IK)  *K*1 *4 ) 

PR  TM  1i5.« 

PRTNI  600 

PRIM  ?O0,  (<|JMWP72  (J)  »JsJ  .4)  »  { A VGtoh* /ii  C .|»  *J*1 »4) 

PRTM  T40n 
PH  I |v  1  400 

PRIM  20  0,  (  (ALP172  ( 1  V-)  i J*1  «4)  ♦  <  ALH2  /2  (!•«)*  K«i»<t)fIsl*  45) 

PRTM  -ifpo 
PR  In  l  600 

PRINT  200*  (CM*  172  ( J)  »  J*l»4)  *  (CMA^/2(H)  *K«1  *4) 

PRTM  3Poo 
PRIM  600 

PRINT  2nO,<«lJMAi72<J)»J*l*4)*  (SUMA2/2  <K>  *K»1  *4T 
PRIM  4000 
PRTM  600 

PRIM  I  20  0* ( A  VGA 1 72 (J) » J*1 »4>  » ( AVGA272 C  K  > »K*1 *4 ) 

PRINT  4^00 
PRINT  600 

PRINT  ?P0,  («CFva13  ( J)  » J*1  »4)  f  (S0ltVA^3tN)  ♦  K* 1 »  4 ) 

PRTM  4400 
Pp  TNT  6o0 

PRINT  ?00*  UHE'/«17?(If  J)  »Jal»4>  ♦  (DfcVA272  (  T  »*)  *K*it4)  •  1*1  »45> 
PRINT  /+f0O 
PRTmT  600 

PRTM  ?00«  ( fl  VGA  i  3  ( «i )  fj=l«4)  »  (AVGA23(N)  .n*1»4) 

PRINT  4°oo 
PR  TNT  6^0 

PRINT  ?O0*  TRG172C*) , <HG272<M »N*i»4) 

200  FORMAT (RFie.f) 

300  FORMAT  (T  HI  ,29X,*U24*t5<>X,*!J24STAR«) 

400  FORMAT  (T  HI  f  ?RX,«V24«*t56Xf*V24STAH*> 

500  FORMAT  (T Hi  ,;?9x.*W24*»5t>X**W24STAK*> 

600  FORMAT  (2  (9*.«3HM8«t J  IX  »«65M8<N  1  nX .* i  i Oni»**  10X **1  9«Mb«  >  */> 

7fl0  FORMAT (TH0t?7X**SLMU24*»55X.*AVGU24«) 

800  F0RMAt(1H  ,//,2 »X»«SUMV24**5bXt*AVGv244) 

900  F0°MAT(TH  ,//»?7X«*SUM2a«,55X»«AVgw24»> 

950  FOHMAT  ( 1  H  ,//»?0X»*SUMWP24».55X*«AVoI»P24*) 

1200  FORMArpHi  ,2?y  ,*ALP1«4**55X»«ALP224*> 

1 JOO  FORMAT  »lHf)t23Xt*  POINTS  Taken  <mI2a>*»37X,«  POINTS  TAKtN 

XA??4)*) 


A-21 


TABLE  10  C0NTINUE0 


i4Qo  format  nHf', -any  «*stM  of  ai?4**43X  »*sum  or-  A224«i 

160ft  FORMAT  rlHr,?py.*AVEMCE  UF  Al?4«  *4  ]  a»*avERAGE  OF  m?24*) 

160ft  FORMAT  (1H1  *?7x*«»A1245TAR*«S4X»*A224SIAh*) 

WOft  FORPAT  (lHft,?-jY,*SLM  OF  A124ST7R  SQUAhEn# , 37* » «SIM  OF  A224Sl*R  SUUA 

xrfo*) 

1B00  FCRPAT (iHo»?lX»*AVERAGE  OF  SUM  CF  A l 24? I AR* »36* »*aVER AGE  0*  SUM  OF 

*  A224STAR*) 

l*Oft  FooMAt  oHo,?7x,*RG12A«»5Sx,*RG??4*) 

200ft  FdRmaT (1H1 t;;9X,*U4rt**56Xf  »U4BS7AR*) 

2i0ft  FORMAT (1  HI ,?9X .*V48**S6X*#V48sTAH*> 

2<;0ft  FORMAT  ( a  HI  ,*W4«**£6X**W4  8s7Ah«) 

260ft  FORMAT  MHP#?7X,*SLM04e*#b5X,*AVGU4H*> 

270ft  FORM  AT ( 1 H  , / / ♦ ? I X , *SUM V4B* , 55X»* AV(,v48* ) 

2»0ft  FORMAt(1H  ,//,?fX»*SUNW4B*»5SX»«AV(iw48#) 

2»6ft  FORMAT  (1H  t//*?bXt*SUMWp4e#«55X»*AV(i»lp46*) 

330ft  FORMAT  (1  Hi  ,?7X,*ALP14e*»S>?X**ALP24b*> 

3b0ft  FORMAT  (lHft,?-aX»*  FOUNTS  TAKEN  ( Al4«>  **37X.*  POINTS  TAKt6  ( 

XA24fi)*1 

370ft  FORMAT (iHof3fX,*SLM  OF  A  14fi« *43x * *SUK  OF  A?4B«) 

390ft  FORMAT  nHr  ,?pX.*AVERAGE  OF  A140#»41A»*ev£RAGE  OF  «24e*) 

4j0ft  FORMAT  (1  HI  ,?7X,*A14@STAR*,F4X,«A24B!sIA>*) 

410ft  FORMAT (IHc ,?3X,*SLM  OF  A14RSTAR  SGuahED* » 37* oF  A248SIAR  S^UA 
*REn*) 

4b0ft  FORMAT  X,*AVEHAGE  (JF  SUM  OF  A146SI  AR**16X,*avERAGE  0m  SUM  OF 

X  A?4HSTaP*1 

470ft  FO°M AT  (lHr  *57x»*RC-148«f5!?x»«H6248*> 

23 Op  FORMAT (1  HI ,?RX»*U7  2*»56Xt*u72s7AH*) 

240ft  FORMAT  OH]  ,29x.*v72*,56Xf«v72STAR*1 
260ft  FORMAT (1  HI ,?9X,#W72*»S6X»»W72STAH*) 

290 ft  FORMAT  OHft  t?7X. *SLMU 72 « »S5X »*AVGU72*) 

300ft  FORMAtnH  ,//,?  'X»*SUMV'2*»55X««AVGv  (2*> 

310ft  FORMAT  ( 1  H  ,//,?fXf#SUM*7<!*»55Xt«AVG,«lT2*> 

315ft  FORMAT  ( l  H  »//t?t>X«*SUM»lP72#*5SX*#Avi>i»P72#) 

340ft  F  0Rmat  n  Hi  t?7x**A(.Pl72#f55X»#ALP2^2«) 

360ft  FORMAT (lHftf?3x»«  POINTS  1 AK£N  (a17?) *»37*» #  POINTS  TAKtN  ( 

XA27 ?)*) 

3B0ft  F0Rmmt  oHft,‘*oX**SLM  OF  A  1 72#  *  43  X  t*SuM  OF 

400ft  FORMAT  MHft»5ex.«*AVEHAGE  UF  A 1 72* * 4 l A t *A xERAGE  HF  «272*> 

440ft  FORMAT (THI,?7Xt*A172STAR*»54A,*A272Sl A**) 

420ft  FORMAT(lHn,p3x,*$LM  OF  Al?2STAR  SGUAHEr>**37x ♦  ««UM  OF  A2?2S»AR  SWU* 
XHEO*> 

460ft  FORMAT (1 Hft *21 X ,*AVERAGE  OF  SUM  OF  Ai F2S 1 AR«»36X»**VERAGE  OF  SUM  OF 
X  *P72STAR*1 

4«0ft  FQRmAt (iHft,97Xt*RG  j ?2« f 56X , *HG272* > 

STOP 
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11*16.1 7,COFY«Bf ( TNPU I *CUTPlT) 
11*16.19. OF»  n000?304  WCROS  - 
11*16. 19.MS  3584  WORDS  ( 

1 1 *1 6. 19 .CPA  .113  SEC. 

11.16. 19. io  .239  SEC. 

11*16.19. CM  2.430  KWS. 

1U16.19.SS 


11 *16.19. PP  1.926  SEC. 

11.16.19.EJ  FNp  CF  J08*  ** 
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IILIST28  ////  EnO  OF  LIST  //// 
ULIST213  ////  ENO  OF  lIST  //// 


o  ooo 


TABLE  10. 


CONTINUED 


PROGRAM  FAI.L  (INPUT  , OUTPUT  ) 

OIMEnSKN  U?4( 1 38*4) *L48 (99*4) *U  <2 (4b,4>  *V?4 (1 *P»4) * V48 (99,4) » 

V7?  (  45  *4 )  ,*24  (138*4),  t,48(*9,4)»W72(45«4)  , SUMU24  (4  )  » 
S(MU40(4)  »SLMU/?(4>  ,SUMV2a(4)  ,SUMV48  (4)  »SUHV72<4)  , 
SI.PW?h(4)  ,9LMVi48<4>  *SUMW72<4l  *AVQU24  u>  ♦AVGU48<4>  * 

AVGU72 (4) »AVGV24  (4) ,AVGV4H (4) »AVGV7?  (4 ) »AVC-W24  <4j  * 

A\/GW4«j  (4  )  ,AVGh  12  (4)  ,DEVU24  (  l3o»4)  ,0£Vl'4B  (99 ,4)  * 
0FVU72(45«4)  ,  Otv/ v24  (1 3d  ,4 )  * Op  v V4ft  ( ?9 .4  )  ,D£vV?2 (45 ,4)  » 

OF  VW?*  <  1  38  ,4  )  »L>E  V*4fi  (99,4  )  »DF  v«72  (45  *4  )  »ALP  124  ( lja»4  J  » 

ALP  140(99,4) » ALP  1  72(45,4) ,ALp224 < 138 ,4) *  ALP 248 (99*4) » 
ALP272  (45,4 )  » SUP 4 1 Z.  (4)  »SuMAi48(4)  »SUM«172(4)  »SUm*224U> 
,cuvfl24e (4) ,5L*A?7  2  (4 )  *AVUA1?4 (4) » AVpa i48 (4) ,AVGa  l72t4) , 
AVGA224 (4) »AVGA?48 (4) ,AVGM2?2 (4) ,UEVAl24 ( 138,4) , 

OFy  A  1 48  ( 99*4  )  *UEvAl72  (4*  ,*♦>  ,0C-VA224  (130*4)  * 

OFV A?<*8  ( 99 » 4  )  »UEVA272  (45.**)  *buEVAll  (4)  »80EVA12(4)  * 
5CEVA13U)  *S0EVA?1  (4)  ,SUfc*A2;>(4)  »S0EVA23<4>  *RGl24  <4>  * 

RO 1  4h  ( 4 )  ,HGl72<4)  ,«G224(4)  ,Rg24B(4)  , PC-*  12  ( 4 )  ,CM  a i24  ( 4)  * 
CN  T  A 1 *8 ( 4 )  *CN7  Al  72(4)  ,CnTm224  (4)  .LMA24O  (4)  » irNTA*  72  ( :+>  * 
AVGA| 1 (4) t  AVGA 12 ( 4 ) , AVGA ) J ( 4 ) »AVGA2l f4) *AVQA22(4) * 

AVGA?  3(4) 

DIMENSION  SAVFU?4 (138,4) ♦ S A VEV24 ( 1 J» » 4 ) »SAVfc W?4  (  i 38 ,4 ) ,SU**P24 (4) , 
X  si8W?<*(  138*4)  »AVGWP24«4)  •  Ut V*h24  ( 1  38  * 4  >  * 

X  s*VEu«»8  (99,4)  »SAVEV48(VV,*^»;>AVEW48{99»<»)  *5UM9P4e(4>  » 

X  91  RW40  (99,4)  ,  A VGWP48  ( 4 >  *DtV#P48 (99,4 )  *  5A VEli 72  ( 45 , <* )  » 

X  SAVEV  /2  (4b»4)  ,SAVE#?2  (43,*)  »SuMWP72  (4)  » SIJBW 72  ( 46, 4)  ♦ 

x  AyGWp72 (4 ) ,CEVWp72 <45,4 ) 


CUTOFF  TS  THE  VALLE  (IN  HAOIANS)  FOh  *m1CH  ALPHA  X 
AND  ALPHA?  CANNOT  BE  EQUAL  TO  Ok  GptATEH  ThAix 


DATA  CUT0FF/#34/ 

00  1  Ja  1  »  4 

1  RCAO  J  00, (U?4 (I»J), I*l»l38) 
DO  2  J  =  1 ,4 

2  RC AO  lno, (U4H(I,J) ,Ial*99) 
DO  3  Jsi ,4 

3  RCAO  mo.  (U72(I,J)  ,1*1,45) 

DO  4  Jsi, 4 

4  rcao  lro,  (V?4(I,J)  ,1*1,138) 

DO  5  Jsi , 4 

5  RcaO  lllfl,  <V4»(I  ,J>  ,1*1  ,99) 

DO  F  Jsi, 4 

6  RtAO  mo,  (V72(I,J)  ,1*1,45) 

DO  7  jsi, 4 

7  Rt-AO  loo,  (W?4  ( I  ,J)  ,1*1  ,1  38) 
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DO  8  J*i  ,4 

8  pe-AD  (W4h(  I » J )  «  l*l»99) 

DO  9  jsl » 4 

9  PcaO  lf>0«  <*72<I  *J)  *1*1  *45) 

100  Format (pFl P.3, 

xM I 7  IPskSFFCS  ThE  INITIAL  NUMBER  OF  HOlwTS  TO  FALH  ELEMENT  OF  THE 
CM  ARR«VS  mSSOCIATEU  WITh  ALPhaI  mixO  ALPHA?  FOR  THE 
•  F»RpE  DTFFEnEFT  TFETAS 

CALL  *MTT(_4.n8.«CNTA124) 

CALL  XMTT (-4 , 13H. .CNTA224) 

CALL  xmtT (.4 »9p. ,CNTA148> 

CALL  Xmtt<-4.Rp.*CNTA?48) 

CALL  xMTT(-4f4s.*CNTA172) 

CA|.l  XMTT  (-4.4S..CNTA272) 

15  FLAGleO. 

FLAG?*  n . 

Call  XMTT  (-4  * (i , SUMJ24 ) 

CALL  xmtT(-4,0»sUFV24) 

CALL  xmtT (.4 ,G.sUFW?4) 

CALL  xmtT ( -4 f r f SUM WP?4 ) 

DO  lo  J* 1  *  A 

C  LOOPS  TO  TFRCiioH  1?0  COMPUTE  THE  SUMS.  AVERAGES*  AND  STABS 

NFFrEn  TO  ruMFUTt  THt  REtD-GtRMANS  FOR  THETA»?a  DEGREES 
C 

t-u  To  I* T  .  1  38 

?!’MU?4 ( JT»SLMU24 (U> *U24 (I . J) -SAVEU24 ( 1 »J) 

<?(  MV?4  (J)*SLMX2a  (^)  ♦V24(I*U)-SAVFV24(l*vi) 

?! PW?4 (J) *SLMW24 (g) *W24 (I *J) -SAVE *24 (1 *j) 

SI  MWP?4  < J> s«UMWF24  L  J ) ♦ *24 ( 1 « J > -SuoW?* U » J) 

io  continue 

Oc  ?0  J* 1 . 4 

AVGU94 ( J) *?UMU24 ( J)/  CNTA224 ( J) 

A  VP  V?  A  |  j)  xciimV?4{  J)  /CNT  A  124(d) 

AvPW?4{j)  *si|mW24  (U)/C(XTAl24  (J) 

AVPNP24  <i„  )  aSuM^PP1*  ( /CNT  A  224  (J) 

?o  Continue 

Do  30  Jal*4 

ru  30  1*1. 138 

TF  (PAvFii?**CI»J)  .F.U.  0)  DE  VU24  ( I .  J)  *U24  ( I  *  J)  »AVGU24  (w ) 
TF(SAvEV?*.(I.J)  .EG.  0)  DEVV?4n.U)*V24(T*J)-AVGV24(U) 

IF  (SavFupa <1 .  J)  .Eu.  01  DEVM?#  II *j)*W24 (T  »J)«AVG*24 (Jt 
I F ( SUPH?4 ( I » J )  .EC.  0)  0EVWP?4(I,J)*«24(T*J)-AVGWP2*(U) 

30  C0NT1NUF 
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OC  40  J*1 *4 

LOOPS  40  ANn  50  COMPUTE  ALPHA1  ANU  alHH/»2»  OMOpe>Ir*b  THOSE 
points  gpeaTer  than  cutoff 

00  40  1*1  *  1  18 

IF ( ALH 1 ?4 ( J  » J)  .EC •  3.15)  GO  TO  40 

AI.PJP4  (I.  J>  *ATAN  <CEVW24(I* J) /UfcVv24  (I*  JM 
Tf(AB«(Alh124(I,J))  .LE.  CU[ohM  b 0  Tu  4r. 
SAVEW?«(I*J)-*424(T«J) 

S  Av  E  V?H  (ItJ):*V24(I*J) 

ALP1?4(I»J>*3.15 

CNT A 1  pH  <i*  )  *CNTA124  <  J>  -1 . 

Fl AGl *i • 

AO  CONTIMJF 

DO  50  J» 1.4 

DU  So  I »1 . l3b 

IF<ALP??4(T,J)  .EC.  3.15)  SO  To  bO 

Al  P??4  < I, J) * AT  AN (DEVWP24(I»J)/UEvu24(lrw)  ) 
Tc(APc(ALH?24(I,j))  .LE.  CUTuHF)  bO  TO  b?j 
Sllf?to?4  ( I  *  J)  *424  ( I  *  J) 

S  AVpl  l?*»  ( I  « J  T*U24  { I  »J) 

ALF224(I.J)*J,i5 

CN'TA??*»(iv)aHMA224(J)«l. 

FLAG?* 1 . 

50  CONTlMiJF 

IF  ANY  prlNT«  HAVE  BEEN  CHOPPED*  RECuMPuiE 

IF(FUAG1  .FC.  I  .CR.  FIAG2  .EQ.  1)  bU  TU  15 
DO  60  jxl ,4 

DU  60  1*1  *  1  38 

IF ( ALP \?4  ( I . J)  .EG.  3*15)  GO  To  65 

Sl)wAl?4(iv)  *SLMA124<  J)  ♦  ALP  124  (  J  ,  J) 

65  IF  (ALH??4 ( I » J )  .EC.  3 . 1 5 »  GO  TO  60 

SU4A224<Os9LMA224<J>  *ALP224  ( 1  *  J) 

60  CONTINUF 

DO  70  J* 1 « 4 

AVGA127  (I,;)  8SUMA124  {  J)  /CNTA124  (  J) 

AVG6?24  (IW)*SUHA224  (J)/CNTA?24  (J) 

70  CONTINUF 

DO  BO  J* 1  * 4 

DO  Bo  f*l *138 

IF (ALP  1 24  < I » J)  .EC.  3.15)  Go  *0  Mb 
HFVA124U, J)*ALP124(I  J)-4VGAl24(J) 

8b  TP(ALF2?4(I*J)  .EG.  3,l5)  GO  10  By 

OFv  A2?*»  ( I  ♦  J)  *  ALP?24  ( I  *  J)  -AvGA?24  <  J) 

80  CqNTINUP 
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00  90  Js 1  * 4 

no  q*  t  =  i«i-,h 

qrEVAi  i  (j)ScnEV4U  ( j)  ♦DEV4124U.J)**? 

srEvA?l(J)=$0EVA21<J) ♦dev A 224 ( 1 , j) **2 
90  CONTIMIC 

DC  1 10  J  =  1  ,4 

avoai l ( j) =qntva l l ( j) /Cntai?4 ( J) 
aVP«?}(ji  ssDfeVA  21  ('j )  /CM  A  224  ( J) 

110  C0NTlM)P 

00  1*0  >j«1  ,4 

90iq4(J)aavciAll(J)/AV6Ai24<J)*«2 
RU??4  (  J  )  a  A  VGA?  1  < J) /AVGA?24  < J) **2 
120  C0NTIIMIC 
125  FLAGlso. 

FL46?sn. 

call  xMiT(«.4,?,suvyAa) 

CAi_l  xv  rT  (-4  .  C  .SUM  V48 ) 
call  Xmtt<-4 ,c.su^w4«) 

CALL  XVTT  (-4»r«SUI*4P4g) 

LOOPS  1*50  Tt-POUuH  2ao  COMPUTt  THt  SUMS.  AVERAGF* »  ANC  STAPs 
•'•FFrEO  T C  CUMFUTE  TFt  REEO-GtpMrtNS  »-0R  IHETAs4p  DEGREES 

00  1  JO  J  =  1  ,4 

no  1 30  1*1.99 

SI  MIJ4C  {  J)  =SLMU4g  <  j)  +(148  (  1.  J)  -SAVFU48(  1  «<J  ) 

SI  MV4P  < J) *SLMV48  <W>  *V48 ( I . J> -5AVFV48 ( 1 
SI  MwaA  (  J>  =SLM*l4fi  U)  *W4H  <  I  *  J)  -sAvt«i4»  (  1  »  j ) 
si  P*P/48  (J)aSUMWp-te  ( J)  ♦  w4p  ( I,  Ji  -SUoW48<  1.  J) 

I3O  C09TIMIF 

Oo  1*0  J»1.4 

AVPD4  0 I J) aSUMU4  0 ( J) /CN7A248  (  J> 

AVP'/4B(  J)  rqUMV48  <  J>  /CNT  A  14B  (  J) 

AVPW48  (  J)  sSi|MW48  ( J  )  /  CNT  Ai48  <  J) 

AvPWpAr  )  *SUMteP*B ( J) /CNTA?48  ( J) 

140  CONTIMIF 

DO  130  J*1 .4 

00  lrf>  Tal.qy 

Tp  (SftVEU48(I.J)  ,EU.  0)  DEVU*eU,J)slj48{T,J>-AVGU48<vi> 
IF (sAvEV4H ( I . J)  .£u.  0)  0EVV4a l I . J) *V4« ( T . J> -AVGV4B ( j) 
TF  (saVE^AH  (I  .  J)  iEU.  0)  DE'''«4b(I»J)sW48(T.JJ-AVG'»48(w) 
TF  (SIIBW4M  ( I  ♦  J>  .EC.  0)  OEVWP48  1 1 ,  J)  «W4«  (  T  .  J) -AVGWP48  (J) 
1*0  CONTIMJF 

DC  1 60  j*l  .4 

LOOPS  1 A ()  AKC  170  COMPUTt  ALPHA  1  ANu  A|kHA2,  CRCmmING  THOSt 
POTxTs  CPF  a ] Ep  THAN  CUTOFF 
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DO  If  0  T  a  1  .9* 

IF  ( ALM  48  ( I « J )  .EG.  3.15)  0  0  TC  160 

Al  Pl4*(T,J)sATAN ( D  EV  w  4H ( I . j) / ut  V  v  4  8 ( I t  J ) ) 
If(4p«(AI.p]48(I»J)),LE.  CUTQf-H  Gu  To  It* 
SAVEW4«  (  I  ,  J) =448  (  I  ,J) 

?AVEV48(I. J)=V4«(I*J) 

ALF148<ItJ)a3.1b 
CNtMaMO  =  CN  Ul48<  J)-l  . 

n  ac-isi. 

1&0  CONT I M  IF 

DO  170  ■  ja 1 ,4 

DO  170  I* )  .9* 

IF (ALp?48(I.J)  .EG.  3.15)  Gq  Tq  170 

41  P?4fl  < I. J) SATAN (CEVWP48  < I  .  J)  /L)EVU48 ( 1 . J) ) 

if ( ar« <A(  (-248 1 1 , j)  >  .le.  Cutoff)  gu  to  nr 

S(JP*48 ( I «  J ) =*48 (I , J) 
SAvFlJ4tt(I.J)=U48(I*J, 

ALP24H(I.J)=3.lb 
CntA?4‘>(iv)=CNTA24«(J)-1  . 

FLflC-2*l  . 

170  CQNT1MIF 

TF  ANY  DO  I NT «  HAVE  BEEN  choppeo*  hecuwpdie 

IF1FUAGT  , FC ,  ^  .C«.  FLAG2  .E<3.  1)  uU  Tu  l?b 
00  1«0  j=l,4 

no  1*0  Tal.99 

IF ( ALP Y  48  t!«j)  .EG.  3.15)  GO  TO  18b 

«?l  HA14*  (  1)  =SUMA148  <  J)  *ALP148  U  t  J) 

1»5  IF  (ALHD48  n.j)  .EG.  3.15)  GO  TO  1«0 

SI.MA24P  ( J)  scuMA24e  ( J )  ♦ALP248  IliJ) 

)H0  CONTIMJF 

DO  1 *0  j*l ,4 

AvDAl 4P pj) =SUPA I48 ( J) /CNTAI48 ( J) 

AVGA?4B  <iw ) aSUMA248 ( J) /CNTAP48 1 J) 

1*0  CONTIMJF 

DO  ?AO  >J*1  ,4 

210  I a ]  *  P * 

TF  <  A|_p  ] 48  <  1 1  J)  .EG.  3.15)  GO  10  ?i5 
DFV/AUtt ( I t J) a ALP  148 ( I ♦ J ) - A vG A  1 48 (J ) 

21^  !F(ALP248,ltj)  .EG.  2.15)  00  |U  ?i 0 

DF^APAb ( 1 ♦ J) sALH?40 (It J)-AVGA248 ( J) 

210  CONTIMJF 

DO  2«0  j*T,4 

ru  92o  t*i .9* 

SrEVAl? ( JJ *5 DEV A 12 ( j) ♦DE*A14M( 1 1 J)**2 
«FEVa?2  {  J)«SOEVA22  <  J)  ♦DEVA248  1 1  t.D*#2 
2<?0  CONTIMJF 
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230 


240 

2t»5 


DC  2-»0  J=1  .4 
AVPA12 ( j) 
AVGA?2  < J) 
CONTIMJF 
00  2*0  >J«1 ,4 
R61 4p ( J) = 
RG24A  ( J)  = 
C0NT1MJF 
FLAGl*^. 
FL4G2»*. 

CALL  XmTT<-4 
CALL  XmtT(-« 
CALL  XmTt(-4 
CALL  XVTT(.4 


=  «;L)tVA12  0w)/CNTAl48  (J) 
s sot V A ?2  M /CNTA248  < J) 


A VGA 12 ( J) /AVGAi 48 ( J) ««2 
AVGA22 (J)/AVGA248(J)*»2 


.C.SUMU72) 

»0.SU*V72) 

*  0  •  SUn-  W72 ) 
,C,SUMWP72> 


LOOPS  3*0  THROUGH  360  COMPUTE.  THE  SUMS,  AVEMAGFS.  AND  STARS 
ixFfrEP  TO  CUMF ijTE  THt  REfcO-GtRMANS  FOR  |HEtA*Y2  DEGREES 


DC  250  J® 1 ? 4 

DO  ?e 0  T*  1  »4b 

Sl’MUT?  ( J)  =SLMU72  ( U )  *U7  2  { I  ,  J ) -SA  VEu7?  ( I  »w  > 

SLMV75  (j)=SLmV72  (j)  ♦V?2U»J)*SaVf  v?2(I*.j) 

Si  MW7?  ( J)  »SIMW72  (0)  *W72<  I.  J)-5AVn»72f  1  *  J  > 

Si MWP7? { J) =SUMWP  ?2  < J) ♦w72 ( I , J) -Suow72 ( I, j) 

250  CONTIMJP 

DC  26A  'J«l  » 4 

DO  p*0  T a  1  * 4b 

AVGU7? ( J ) sSL  MU72 < V ) /CNTA272 < J) 

AVGV7?  (  J)  *SLMV72  ( >j )  /CNT  A 1  72  ( J) 

AVGW72  ( J)  =SLM»i72  (sj)  /CNTAl7£  (U> 

AVGi-'P7c  (  J)*SUMWP72(  J)  /CNTA2(2«vl) 

260  CONTINUE 

DC  270  j*l  ,4 

CO  27 0  T*1 *46 

rP(SAVEU7<r(I,JJ  ,EU.  0)  DEVU72(I,U)all?2(T«J»-AVGU72U) 
r*(SAvEV72(  I*U)  .EG.  0)  DEVV7*<I,J)3V72<T.J)-AVGV72<,j) 
IP  (SAVEW72CI»J)  .Eu.  0)  DEVW72U»U)=W72(T»U>-AVGM72(w; 
TP  <SUPW7?(I»J>  .EC.  0)  DEVWP72U’,u)=W72(T,J)-AVGWP721J) 

27o  continue 

LOOPS  2P0  AND  ?S*0  COMPUTt  ALPH«1  Anu  Al.PHA?.  CROPPING  THOSt 
POINT®  GRFATfcH  Than  LLiIOFF 
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00  2B0  js] ,4 

DO  bPO  T  a  1 ,4b 

IF(ALMT2<T,J)  .EC,  3. IS)  00  TO  26o 

Al  Pi  75  < I , J) =  AT  AN VQEVW72 (I, j)/utVV/2(I,j) j 
IP  (AR?  (AL^172  (I,  J)  )  .LE.  CurOf-h)  GO  TO  2 8o 
S«VEW72CI.J)-W'2<I.J) 

SAVEV72(I.J)sV72<I*J) 

ALP  17? ( I ,  J) =3 . lb 
CN'TA17?('v)aON  IA1T?(J).1, 

FL«(?1  =  1. 

2Bq  CON'TLMJP 

DO  2*0  j=1 ,4 

PO  ?9 0  Ts  1 » 4b 

IF  ( ALp??2 ( T , J )  .EG.  3. IS)  Go  To  290 

At  P?7?(I,J)a&lAN  (CEVWP7?(I,J)/1)Evo7?(I.w)  ) 
IP  <AP«  (ALP2?2d»J)  )  ,LE.  CUTOFC)  GO  TO  290 
SIJPW7?  (I ,  J)  =h?2  <  I ,  J) 
CA\/EU74:(1»J)sU72(I,J) 

ALP27?(I»J)=3.i5 

OtA?72(i,.)sCM  A272<J)-1. 

PLSOPal . 

290  CONT IMJP 

IF  ANY  PPlNTc  HAVE  5EEN  CHOPPED,  HECuPPulE 

IF(FLAG1  .PC.  1  .CP.  FLAG2  .EO.  1)  GO  TO  255 
00  j=l  ,4 

DO  3l 0  I a  1 *4b 

IF ( ALp I  72 ( T , J )  .EC.  3.15)  GO  TC  3lb 

51  HAi 7? ( J) =SUPA172 ( J) +AIP] 72 ( i * J> 

315  IF  IALH?;2(T,j)  .EC.  3.15)  GO  TC  3 1  (i 

pt  MA?7?  < J) r PUMA 2 72 ( J) ♦ALp2 '2( 1  ♦  J) 

310  CONT1MJF 

DO  3«n  j  =  l  ,4 

A VGA  172  Cv) SSUMA 172 ( J) /CNTAl 72 ( J) 
AVga?7?(iv)*guMA272(J)/CNTA272(J) 

320  CONTIMJF 

DO  3J0  ij*1  ,4 

DO  330  Ts 1  ,4b 

IF  ULF17?U,J>  .EG  »3 • lb)  GO  lu  33b 

HFVA1 72  (1,0) s ALP  1 72  1 1 , J ) -« vM , 72 ( J) 

3^5  TP (ALC2721I ,J)  .EG*  3.1s)  00  10  330 

PE'< A ?7?  ( I ,  J)  «ALP272  ( I , J) -avGA? 72 ( J> 

330  CONTlMjP 

DO  3*0  J=7 , 4 

DO  3*0  Ta]  ,4t> 

prEVAi 3 ( J) SPQEVA 13 ( J) ♦OEVA 1 72 ( i » J) **2 
FrEVA"3 { J) s« DEV A 23 ( J) tUEV 4272 (I, J)**? 

3*0  CONT  I M  IF 
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00  3»0  j»l,4 

Ai/GAi 3 ( *RUtV A  1 3<w>/LKTAi72(J) 

A VP  49 3  (  _n  *SDtVA?3  (w  )  /CNTA272 ( J> 

360  COMIMJP 

00  360  j*l  .4 

RG17?  (J)aAVfiAl-a(J)/AVGAl72(J)**? 

OG???  <  J)  sAVr,A23  (  J) /AVGA?72 ( J) *«2 
360  C0NT1MIP 
PPIM  ISO 
PPIM  AO  0 

PPINl  200*  ( (D?4 < I , J) » J® 1 ,4) » (DEVU24 (I ,6) ,K*1 »4) » 1*1 ♦ 138) 

PPIM  400 
PRIM  600 

PRIM  200,  <  <V?4(I,J)  »J»1,4)  »  (DEVV2?*(l,iV)*K*l»4)»l*l»l38) 

PRIM  ^00 
PRINT  600 

PRINT  200,  (  (W?4  (I  ,  J)  ,  J*1 ,4)  *  (0EVW24  ( I  ,N)  ,K«tl  *4)  *  1*1  »  138) 

PRIM  700 
PRtnT  6*0 

PRIM  900,  (SUMU24  ( J)  »  J*l*4)  »  (AVGU?**  (K,  ♦K*l*4) 

PRIM  o0n 
PRIM  600 

PRINT  ?00,  'SUMV24  ( J)  ,  J*  1,4)  ,  ( AV6V2MK)  ,K*1  ,4) 

PHIM  ROO 
PRIM  firp 

PRINT  200,  (SUM*24  <  J)  ,vl=l  ,4)  »  (AVGW2*  U)  ,K*1 ,4) 

PRIM  Rc0 
PRT^T  Ao 0 

PRIM  ?0  0,  («UM«‘H24  (  J)  *  J® 1  *  4)  #  (AVGi KP*4  (J) »J*1»4> 

PRUT  1300 
PRTNT  600 

PRINT  200,  (  ( AL° l 24  ( I  »iw  )  »  U®  1  »4)  ,  ( ALPC24  <i,K)»K*l,4),I*l»138) 
PRIM  1-sOO 
PRlNl  60  0 

PRIM  200,  (CNTA124 (J) »  J® l ,4 )  t (CMA224 (k) »K*1»4) 

PRIM  I4C0 
PRIM  60  0 

PRIM  200, ( eUMA l 24 <J) »J«1»4> » ( 5UMA224 (k) »K*1 ,4 ) 

PRTM  1®00 
PPINI  6p0 

PR  INI  200, (AVGA124 (J) »  J® 1 »  4 ) t ( AVGA224 (K ) ,K*1 » 4  > 

PRIM  l^oo 
PRIM  600 

PRIM  ?00, («CEVAll (J) » J»l*4>  * (S0EVA21 (K) fK*I»4) 

PRTNT  1600 
PRlNl  600 

PRTNT  900,  (  (DFVA124  ( I , J) f J*1 *4) ♦ (Dt vA2?4 ( I ,6 > *«*l *4 ) * I«1 • 138 
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PRIM  1  0  0  n 
PRT*I  400 

PR TM  200.  (/JVC-All  C,J)  ♦'wal»4,  .  (AVGA2)  (K)  ,4) 

PRINT  icpo 
PRINT  600 

PRIM  20  0,  (is. )  *  J=l«4)  ,  (Hf5224  (N)  ,  Kb i *  4 ) 

PRINT  2000 
PRINT  600 

PRIM  200  ♦  <(U4MU,  J)  *  J*1  *4)  *  (OEVU4«(I,M  ,K*1  *4)  »  lsl,99> 
PRINT  2’ 00 
PRINT  600 

PRIM  200*  (  <V4H(I»J)  *Jsl*4)  *  (DEW*a(T,M  ,K*1*4)  *1=1*99) 
PRINT  ?,0f 
PRTnI  600 

PRTfjl  ?flO.  (  (W*b(I*J)  *J»1*4)  *  (0EVWaoU*A)  *K«I*4)  *l*i*99) 
PRIM  2*00 
PRTNI  400 

PRINT  200,  (SUMU4P { J) » jal *4) » (AVGUAB (K> »K*1 *4) 

PRTNI  2700 
PRTNI  *r0 

PRINT  200*  (?l)MV4e  (  J)  *  Ja  1,41  »  (  AVGV4BIK)  *K  =  1  *4) 


♦U.S,  GOVERNMENT  PRINTING  OFIICE  :  1980  0-721-894/288 
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